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Nonperturbative QCD approach is systematically derived starting from the QCD La- 
grangian. Treating spin effects as a perturbation, one obtains the universal effective 
Hamiltonian describing mesons, hybrids and glueballs. 

Constituent mass of quark and gluon is calculated via string tension. The resulting spec- 
trum of mesons, hybrids and glueballs obtained is in good overall agreement with lattice 
(f) ■ data and experiment. 
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1 Introduction 

Hadron physics is a well-established field with enormous experimental treasury of facts; on the 
theoretical side there are different models, some of these are well developed and have an im- 
pressive predicting power. The theoretical situation cannot be however considered satisfactory, 
taking into account that it is now clear that QCD is the real theory of strong interactions and 
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can be used for selfconsistent construction of the hadron theory. In addition to that new power- 
ful methods in lattice simulations of QCD yield detailed information not only about spectrum 
of hadrons, but now also about scattering and hadronic matrix elements. 

The reason for this discrepancy was lack of analytic nonperturbative (NP) methods in QCD, 
which could be used for systematic derivation of QCD-based models. Recently the situation 
started to improve, and new powerful analytic approaches have been suggested and checked 
versus experiment and lattice data. Remarkably many of the previous model results are derived 
now by these approaches in a rather direct way from the first principles, i.e. from the QCD 
Lagrangian, and it is gratifying that these features are supported by experiment and lattice 
data. In addition new NP methods allow to resolve old problems, which plauged existing 
models. The (incomplete) list of problems includes: 

1) The problem of Regge slope. Relativistic quark models (RQM) predict for it the value 
(8a)" 1 , where a is the string tension, while the string models give (27T0") -1 , which seems more 
realistic from physical point of view. 

2) The problem of Regge intercept. Both RQM and string models predict too high values 
for the hadron masses, which is mended by introduction of large negative constant in the 
Hamiltonian, of unknown dynamical origin. 

3) Perhaps a more fundamental problem is the problem of constituent mass of quarks and 
gluons. The latter are introduced as a phenomenological input to reproduce spectrum and static 
characteristics and are different in the relativistic and nonrelativistic versions of the model. The 
dynamical origin of constituent masses is unknown. 

4) The problem of quark radius and formfactors. In RQM ( to take a most developed model) 
radius of hadrons comes out to be too small which is usually mended by introduction of some 
effective quark radius, the origin of this radius is not understood. 

5) A very fundamental problem which lies at the basis of the whole QCD, is the connection 
between the quark models of hadrons (based on the minimal number of degrees of freedom 
- those of valence quarks and gluons) and high-energy scattering, including DIS, where bare 
(current) quark degrees of freedom enter, including quark sea and gluon d.o.f. As an example, 
about one half of energy-momentum is transported by a gluon cloud around proton, while in 
the constituent quark model this gluon cloud is not evident. 

In this review we attempt to establish a basis for answering these questions, using new NP 
approaches in QCD. These include: 1) a new background field theory, where NP fields play 
role of the background, and the separation of perturbative and NP is made formally exact by 
virtue of the 'tHooft identity; 2) Relativistic path integral representation for quark and gluon 
Green's functions; 3) Method of relativistic Hamiltonians with auxiliary functions; 4) A very 
recent method of effective contour-dependent quark Lagrangians, yielding a simple picture of 
hadrons in the large N c limit. 

The space limits of this review leave out many details of the formalism, which can be found 
in the original literature. Moreover a very important part of the hadron physics - the strong 
hadron decays - is left over completely for future. 
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2 The QCD vacuum, confinement and chiral symmetry 
breaking (CSB) 

QCD is believed to be the theory of strong interaction. The evidence for that statement 
comes from comparison of perturbative QCD and the so-called QCD-motivated models with 
experiment and also from lattice data. 

There is some gap between models, treating nonperturbative (NP) features of QCD, and 
QCD Lagrangian, which is now filled in by new NP methods and Monte-Carlo calculations. 
It is the purpose of the present review to describe these new NP methods, confronting their 
results with experiment and lattice data, thus demonstrating that one can indeed derive the- 
ory of hadrons from the first principles of QCD with few assumptions, which can be checked 
independently. 

To describe properties of the QCD vacuum and NP dynamics of quarks and gluons one must 
consider the Euclidean QCD Lagrangian (since NP configurations turn out to be Euclidean) 
and the QCD partition function has the following form. 

Z = J DAD^D^ + exp[L + L x + L int ] (1) 

where we are using Euclidean metric and define 

Lo = -\ J d*x(Fp 2 , (2) 

Li = -iJ f i/j + (x)0 + m f ) f ip(x)d 4 x, (3) 
f iP + (x)gA(x) f ^{x)d A x. (4) 



-'int 



Here and in what follows ^ att denotes quark operator with flavour /, color a and bispinor 
index a. 

The quantum theory of QCD is believed to be fixed by equations (1-4) and one prescribed 
scale, e.g. k QCD . 

There are two main features of NP QCD, which make our world and ourselves as we are: 
confinement and CSB. 

Confinement is a property of gluonic part of Lagrangian (2) (for a review see [1]) and most 
easily defined for static quark Q and antiquark Q neglecting light quarks, in which case the 
celebrated order parameter is the Wilson loop (to be often used below), 

(W(C)) = J DAexpL ■ Pexpig J A^. (5) 

We suppress here the gauge fixing and ghost terms for simplicity. For large contours C one has 
Wilson criterium of confinement [1] 

(W(C)) = exp{-aS - 7 P) (6) 

where S is the minimal area inside contour C and a is the string tension, while P is the 
perimeter of the loop C. Note that theory must connect a and Aqcd ~ the only scale of QCD, 
how it is done on the lattice and in some models - see in [2]. 
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Confinement in the form of area law (6) produces linear potential for heavy quark and anti- 
quark at large distances, V^q(R) = aR, and in section 5 we shall find also other contributions 
to Vqq(R) due to perturbative gluon exchanges and spin-dependent parts. It is important 
that light-quark loop contributions are suppressed by powers of 1/N C and can be accounted for 
perturbatively. 

Whereas confinement enters directly via W(C) in the heavy quark Green's function (see 
section 5), it is much more involved problem for light quarks, especially for quarks in pions, 
since their Green's function cannot be expressed readily through W(C). Physically this fact 
is connected to the Zitterbewegung of light quark which is connected to the phenomenon of 
chiral symmetry breaking. To proceed one remarks that massless quarks can be split into left 
and right chiral parts, ty(x) = ^/r + ^> l and ^ l enter Lagrangians (3) and (4) in the 
additive way, exemplifying the chiral symmetry of the group U^Nf) ® Un(Nf). It was shown 
in [3], that nonperturbative gluonic fields responsible for confinement, also produce CSB for 
light (massless) quarks. Concretely the appearance of the string (which transforms as Lorentz 
scalar) for a light quark is already a sign of spontaneous CSB, and the string appears as a 
solution of the nonlinear equation for the light-quark Green's function [3]. 

For most mesons, however, not pions and kaons, chiral and spin effects can be treated as 
perturbation of the basic confinement dynamics. The corresponding Hamiltonian [4,5] looks 
much simpler, and reduces in some limiting cases to the familiar relativistic potential models 
[6,7,8]. We derive these Hamiltonians in section 5 in the cm. system. 

Of special interest in modern hadron spectroscopy are glueballs and hybrids. Here relativis- 
tic hamiltonians can be easily generalized to include gluon degrees of freedom, if gluon spin 
interaction is considered as perturbation. The spectrum obtained in this way is discussed in 
section 6 and is in good agreement with lattice data. One should note, that nowhere in our 
formalism we introduce constituent quark or gluon mass, and the only input mass parameters 
are string tension a, and Aq C d ( in addition to current quark masses rrif, renormalized at the 
scale lGeV). With such input one obtains a wide variety of hadron states and its properties 
in good overall agreement with experiment and lattice data. We also define and calculate the 
constituent mass of quark and gluon in terms of a. 

3 Perturbative and nonperturbative field configurations 

One of basic mysteries of QCD is the fact, that gluon field plays two different roles: 

a) gluons are propagating, and at small distances this process can be described perturba- 
tively, leading in particular to color Coulomb interaction between quarks (antiquarks); 

b) gluons form a kind of condensate, which serves as a background for the propagating 
perturbative gluons and quarks. This background is Euclidean and ensures phenomena of 
confinement and CSB. 

Correspondingly we shall separate the total gluonic field into perturbative part and 
nonperturbative (NP) background B^. 

A^ = B^ + (7) 

There are many questions about this separation, which may be answered now only partially, 
e.g. what exactly is the criterion of separation. Possible answer is that perturbative fields a M get 
their dimension from distance (momentum), and therefore all correlators of fields a M (in absence 
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of Bfj) are singular and made of inverse powers of (x — y) and logarithms, where enters the only 
dimensional parameter of perturbative QCD - Aq C d- Therefore evidently any dimensionful 
constant, like hadronic masses or string tension cannot be obtained as a perturbation series. 
In contrast to that, NP fields have mass dimension due to the violation of scale invariance 
which is intrinsically present in the gluodynamics Lagrangian. The origin of separation (7) 
is clearly seen in the solutions of nonlinear equations for field correlators [9]: a perturbative 
solution of those leads to singular power-like field correlator, whereas at large distances there is 
a selfconsistent solution of the equations, decaying exponentially with distance with arbitrary 
mass scale, since equations in [9] are scale-invariant. Full solution including intermediate 
distances produces mixed perturbative-nonperturbative terms, containing both inverse powers 
of distance and exponentials. For these terms criterion of separation fails. 

One can avoid formally the question of separation principle (and of double counting) using 
t'Hooft identity [10], which allows to integrate in (1) independently over and a^. 



Here L(A) = L + Li+L int , and is taken to be B^+a^. For the exact formalism starting from 
(8) we refer the reader to [11,12,2], and here we only quote the form of background propagator 
Gfj, u of the gluon, which is found from L[B + a) and will be used in what follows, 



where D x = d\ — igB x , and F^ u = d^B u — d u B lM — ig[B^, B v }. Neglecting NP background one 
calculates the perturbative static potential defined gauge-invariantly as 




(8) 



= {D\5^ + 2igF^ u )- 1 



(9) 



V(r) = —Um lnW(r,T) 



(10) 



Recent two-loop result in configuration space is [13] 




r 



(11) 



and in momentum space one has [13] 



V(q) 



-c F 



4:%av(q) 




with 



a v (q) = aMs(q)(l + ai^MsiQ) + ^hs) 
and for the potential in position space one has [13,14] 

a R (l/r) = aMs(l/r)[l + c^Msi 1 / 1 ") + c 2alis) 
where in quenched QCD (n/ = 0) one has 



(13) 



(14) 



ci = 1.832, c 2 = 1.758 



(15) 



Employing the lattice measurement of a(l/r) in one obtains [13,14] 

Ag } 5 = 0.602/ro, r w 0.5/m. 



(16) 
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With this value of Ar one obtains the renormalized or(1/7") which diverges at the Landau 
ghost pole, situated at r = 0.2/m. 

Thus the perturbative potential is unacceptable in this form at r > 0.2/m. However analysis 
of lattice data in [14,15] reveals, that Landau ghost pole is absent in the data. Potential models 
usually postulate saturating behaviour of olr at large r of the type [7] 

(1/r)= U«(one-loop), r<r 
y const = aR{r ), r > r 

In this way one introduces a new parameter aR(max) = OR(ro), usually chosen around 0.4. 

Recently the physical mechanism producing freezing of cxr at large distances was identi- 
fied in [11,12]. It is connected to the fact, that in the confining background gluons cannot 
propagate very far from the sources, and typically this distance is invesely proportional to the 
corresponding hybrid excitation energy (which is of the order of 1.0 -j- 1.2GeV). 

At large N c and in the Euclidean region, q 2 > 0, one can predict this saturation (freezing) 
of a s theoretically to one loop as [12] 

a (q) = — (18) 

where ip is the Euler function, ip(z) = T'(z)/T(z), m 2 is the scale of radial excited states at 
large n 

Ml = nm 2 + M 2 (19) 

and Mq the lowest mass state in the channel. 

When the argument of ip function is large, one has asymptotic representation 

Mz) =lnz- — -J2 -=^k (20) 

Correspondingly a(q) asymptotically assumes the form 

« M (g) = -^, b = l -±N c -\n f (21) 

A naive "explanation" of (21) is that gluon acquires the mass M = m g which eliminates 
Landau ghost pole. 

In reality gluon does not have the mass (this would violate gauge invariance and render the 
theory nonrenormalizable), but due to confinement is connected by the string to the quarks 
and this creates the mass of the whole system - the hybrid. 

The two-loop asymptotic form of a in the position space looks like [2,12] 

/, ; \ 47T - bi In In a(r)-, ■, 

= , i , M + r j ( \ } 22 

b m a(r) b In a(r ) 

where a(r) = M °1~J — , and b x = 102. 

Note that the combination Mq + q 2 , Mq + 1/r 2 in a is the renormalization group (RG), 
invariant, since both external momenta and background fields (gB^), defining Mq are RG 
invariant [2,12] (while g and separately are not). 
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There are many phenomenological arguments in favour of freezing of a s at large distances, 
for a rewiew see [16] and more recent discussion in [12]. A different approach to the freezing 
behaviour, based on some assumed analytic properties of a has been suggested in [17]. 

Summarizing phenomenological situation one can state that all known facts are supporting 
saturating a in Eucledian region at the level a(max) < 0.8. 

Recent lattice data allow a more direct computation of a{l/r) and thereby the freezing 
phenomenon. In [18] the form (22) was used for ap, entering the force between static quarks 
(which is easier to measure on the lattice) 



"IT) - (23) 



The lattice measurements of [15] are compared in Fig.l of [18] to the theoretical curve (23) with 
M = lGeV and = 280MeV. One can see there a good agreement, whereas the nonfreezing 
behaviour (M = 0) clearly contradicts data at r > 0.4/m. 

As one can see from [18], lattice calculations give a clear direct evidence for the freezing of 
the coupling constant, with the maximal value 

a Lattlce {max) < 0.5 (24) 

A recent experimental evidence for freezing of a s was found in the detailed analysis of charmo- 
nium spectra in [19]. 

At this point it should be noted that freezing (saturation) is a property of Euclidean (space- 
like) interactions. In the timelike region (q 2 < 0) Equation (12) yields complex logarithm, which 
is drastically important for time-like formfactors and Drell-Yan processes. From physical 
point of view logarithmic branch points of a s are not relevant, and the analytic structure 
of a s is better displayed by the Equation (18) valid in the large N c limit. In this limit all 
singularities of physical amplitudes are poles, and in the physical scheme of a s renormalization, 
when a s is directly connected to some physical amplitude, a s should have only poles, which are 
clearly displayed in (18). At large \q 2 \,q 2 m 2 , one can average over many poles and return 
to the complex logarithmic dependence in (12), in approximate agreement with standard a s 
expression, where Mq = 0. 

So far we explored purely perturbative expressions for static potential or modified due to 
the NP background. Now we turn to the purely NP static interaction, i.e. to the confining 
potential. To this end one rewrites the NP background part of the Wilson loop using the 
nonabelian Stokes theorem [20] as 

(W(B)) = r^tr(P exp ig ^B^dxJ B = 

= ^-tr(exp ig f d<j^F^(u,x )) B . (25) 
Here S is a surface bounded by contour C and x - an arbitrary point on S, we also defined: 

Fp,(u, x ) = 0(x o , u)F^{u)<l>{u, x ) (26) 
and $(x, y) is the parallel transporter along some contour from y to x. 

rx 

<j)(x,y) = Pexp ig B^zjdz^ (27) 
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To evaluate the last integral in (25) one can use the cluster expansion theorem [21] 



(W(B)) = -Lfrexpjr^- [ da{\)... f da(n) « F(l)...F(n) » (28) 

iV c n=l ^" 

Here F(k) = F^u™, x ), da(k) = tfcr^ («<*>). 

The double brackets in (28) denote the cumulant, or connected average, e.g. for n = 2 

« F(1)F(2) »= (F(1)F(2)> - (F(1))(F(2)) (29) 

For more details and definitions of cluster expansion see [21]. 

Throughout this review we shall consider cumulants or field correlators (FC) ((F(l)...F(n))) 
as the basic NP input, which defines all dynamics of the quark-gluon systems, both for rela- 
tivists or nonrelativistic situations and potential or nonpotential regimes. It is a rigorous and 
explicit language which allows to derive properties of confinement and CSB from the structure 
of cumulants, as will be shown below. On the other hand, FC are gauge-invariant and Lorentz- 
covariant quantities which can be found independently, both analytically and on the lattice. 
Analytic studies of the lowest cumulant (FF) have been done recently in [9], using selfcoupled 
equations which exist for large N c . Lattice studies are by now numerous for (FF) [22-25] and 
also quartic FC-cumulant FFFF > - was studied in [26]. 

Let us discuss first confinement, i.e. the area law (6) in terms of FC. To this end we rewrite 
(28) as 

(W(B)) = l-trexp[~ ^ da pu {u) jf da p<T {u)K„{u, v, C)\ (30) 
where the global correlator A is introduced, 

A pu , P a(u,v,C) = g 2 < F pi/ (u,x )F pa (v,x ) > - 

oo (j a )n , 

- 2 V ^f- / da(3)...da(n)perm < F(u)F(v)F(3)...F(n) > (31) 
„ a J 

and perm denotes sum over permutations of F(u), F(v) and other terms in the cumulant. 
Consider now the lowest term in (31), FF ^>. 

Both from lattice measurements [ 22-25] and analytic study [9] this term decays expo- 
nentially with \u — v\ and the characteristic correlation length T g , which is called the gluon 
correlation length, is rather small, 

T g = 0.2 + 0.3/m (32) 

This quantity is a basic characteristic of the QCD vacuum, defining different dynamical regimes 
for systems of the range R, in cases R(T g or R ^> T g , as will be shown in chapter 5. For an 
earlier discussion of T g see [27]. 

Consider now Wilson loop of large radius R,R^> T g . Then in the integral (30) the 
generic situation is when \u — xo\, \v — xq\ of the order of R. In this case the correlator 
(F tll/ (u,xo)F p(T (v,xo)) does not depend on xo (up to the terms 0{T 2 /R 2 )) and can be rewritten 
as the function of \u — v\. The general decomposition then can written as [28] 

gHr(F lxv (x 1 )^(x 1 x 2 )F Xa (x 2 )^(x2,x 1 )) = iV c [(<5 MA 5 )D{u)+ (33) 
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+-( ~ u x ■ 5 ua + u^5 ua + perm)Di(u)] 




Insertion of (33) in the integral yields readily the Wilson area law 



(W(B)) = exp(-(rS - 7 P) 



(34) 



where a is expressed through the function D(x), while the perimeter term 7P is due to both 
D and D 1: 



The ellipsis in (35) implies possible contribution of higher correlators in (31). For all of them, as 
well as for the Gaussian correlator, a condition necessary to contribute to the string tension a 
is that they should contain "Kronecker component", like D(x), which is a coefficient of product 
of Kronecker 5 symbols. 

As a consequence this term violates Abelian Bianchi identities [28] (which can be checked 
acting on both sides of (33) with operators e af 3^ u -£^) and hence is connected with the contri- 
bution of abelian projected monopoles. This is one of the possible illustration of the physical 
mechanism of confinement induced by correlator D(x) (and similar higher correlators). 

Let us come back to the gluonic correlation length T g , which together with a constitute 
the basic nonperturbative scales. Actually one can calculate T g through a, as it was done in 
[9] using equations for correlators, and through gluon condensate using QCD sum rules [29]. 
In both cases T g <C Aq C d and T g <C which means that T g signifies a new physical scale. 
The physical meaning of T g can be understood, when one calculates the field distribution in 
the QCD string [26], where Gaussian correlator D(x) yields a very good description of lattice 
data and reveals that T g defines the width / of the QCD string (where fields decrease by 50%), 
namely / ~ 2T g . From lattice data [26] one has 



In what follows we shall sometimes use the limit T g — > 0, which we shall call "the string limit". 
Actually this is the limit of thin strings which is as a rule much simpler than dynamics of the 
realistic QCD string of finite width ( the same limit is always assumed in the theory of strings 
and superstrings). 

The finite value of T g brings about specific effects. E.g. for the potential of static quarks 
V(r), which is obtained from the Wilson loop as in (10), one can express V(r) through D(x) 
and higher correlators, and analyze the large and small-r behaviour, which yields [28] 



where Co ~ uT g , c 2 is expressed through D(x), C4 through quartic correlator and so on. Eq.(37) 
demonstrates that NP fields (B^) are soft and yield analytic behaviour of V(r) at r <C T g , 
corresponding to the OPE. Interference of perturbative and NP fields may change this picture 
at small r yielding linear term in r, as was shown in [30]. For more discussion in connection 
with OPE and QCD sum rules see [31]. We shall not go into details of this phenomenon, but 
mention that linear behaviour at small r is phenomenologically necessary for good description 
of fine structure of heavy quarkonia [19,32] and was found recently on the lattice [14]. 




(35) 



Z = 0.4 ^ 0.5/m, T g = 0.2 0.3/m. 



(36) 



V(r) = ar - c , r — > 00; V(r) ~ c 2 r 2 + c 4 r 4 + r <T g 



(37) 
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We end this section with the discussion of the validity of Gaussian approximation, when only 
quadratic correlator (FF) is taken into account. There are arguments in favour of Gaussian 
approximation both from comparison with lattice data and from the structure of the method. 

The first arguments come from the analysis of field distribution inside the string made in [26] , 
as was discussed above. There the largest omitted contribution (from the quartic correlator) 
amounts to less than few percent. Another comparison is for the static potentials in higher 
representations, which demonstrate in numerous lattice data [33] (for earlier data see in [1]) a 
clear Casimir scaling. This is in agreement with the Gaussian approximation, while the quartic 
term would violate the Casimir scaling [1],[34]. From the numerical accuracy of the latter one 
can again deduce, that Gaussian approximation is valid within 10% of accuracy, or even better. 

Finally one can estimate the parameter of cluster expansion as follows. All observables are 
expressed in terms of connected FC of gluonic fields, ((F(l)...F(n))). E.g. the string tension 
for heavy (static) quarks is an infinite sum of FC of the field F u = E 1 integrated over the plane 
14: 

, n n—1 

° ~ II d 2 r i ((E 1 (0)E 1 (r 1 )E 1 (r 1 + r 2 )...E 1 (J2r))) (38) 

One can identify parameter of expansion in the sum (38) to be (only even powers of n enter 
the sum) 

C = {EiT 2 g f (39) 

where T g is the gluonic correlation length in the vacuum defined by the exponential decay of 
FC, and E\ = g 2 {(E") 2 ) = ^G 2 « (0.2GeV) 2 while G 2 is the standard gluon condensate. 

Lattice calculations confirm that T g is rather small [22-25], indeed T g 0.2 0.3 fm and 
therefore ( is a good expansion parameter 

C = 0.04 ^0.1 (40) 

The regime (40) ( C 1 which seems to be characteristic of real QCD, can be called the 
regime of the weak confinement . In this case the dynamics of quarks and gluons is adequately 
described in all known cases by the lowest (Gaussian) correlator. 

One should mention the negative feature of Gaussian approximation. Whereas the global 
correlator A^^ - containing the infinite sum over reduced correlators (cumulants) yields the 
area law for the Wilson law, not depending on the shape of the surface being integrated in 
(31), and hence reducing to the minimal surface,retaining only Gaussian correlator one gets 
as a price for simplification the parasytic dependence on the shape of the surface S in (30). 
Therefore saying about the dominant role of the Gaussian correlator and smallness of discarded 
terms one should specify that estimates refer to the minimal surface. On the other hand for 
an arbitrary surface of any weird shape the contribution of higher correlators can be large and 
this is what exactly needed for compensation of the extra area contribution of the Gaussian 
correlator. 

The situation here is the same as in the QCD perturbation series, which depends on the 
normalization mass fi for any finite number of terms of the series. This unphysical dependence 
is usually treated by fixing ji at some physically reasonable value /io of the order of the inverse 
size of the system. 

In what follows we shall not use Gaussian approximation, expert for calculation of spinde- 
pendent contributions, which amounts to a relatively small correction in good agreement with 
lattice data and experiment. 



10 



4 Background perturbation theory and definition of 
hadron states 

To define the perturbation theory series in ga^ one starts from (8) and rewrites the Lagrangian 
as follows: 

L t ot = L gf + L gh + L(B + a) = 

L + L 1 + L 2 + L int + L gf + L gh (41) 

where Li have the form: 

L 2 (a) = \ u (L>l5^ - D^D U + igF^a^ = 

= L:id?dz% u - D^Df - g rXuK > ( 42 ) 



D? = d x -5 ca + g r ba B b x = b x 



4 

Lint = -\{D,{B)a u - D u (B)a,) a g f abc ay u -\g 2 f abc ayj ae l a^l (43) 

L ghost = -e+iD^D^B + a)) ah 6 h . (44) 
To calculate Green's functions of any hadrons one can use (8) and write 

G h (X;Y) = const J DB^{B)D^D^Da^+{X)^< in {yy Ltot 

= (*l(X)* in (Y)) Bti>>a (45) 

Here ^/,m are initial and final hadron states made of B, if), ijj, a. To calculate integral Da^ one 
can use perturbative expansion in ga^, i.e. neglect in first approximation L int , containing terms 
a 3 , a 4 and take into account only quadratic terms L 2 , L g f, L gh . contribution of L x was studied 
in [11] and it was shown there that for the most processes (e.g. for hadron Green's functions) 
it can be neglected (with accuracy of about better than 10%). 

It is convenient to choose background gauge for a M , D^a^ = 0, and take gauge transformation 
in the form 

B, - U + (B, + l -d,)U, a, - U + a,U (46) 

Now one can choose ^//,^/j n as gauge-invariant forms built from a M , ■?/>, ■?/> and to average 
over 1?^, ■?/>,'?/>, a M as shown in (45). Integrating over DB^ might seem an impossible adventure 
due to unknown rj(B), but we shall see that it can be always written as products of Wilson 
loops with some insertions, which can be treated in two ways: 

i) or in the form of the area law. Then contribution of fields B^ is given and fixed by the 
string tension o~(B). Equating a(B) = a exper one fixes contribution of the background fields. 
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ii) Using cluster expansion, and expressing result through lowest correlators, e.g. (FF), or 
D(x) and Di(x). In this case background is fixed by these functions taken as an input, e.g. 
from lattice data. 

The rest integrals, over Da u DipDip are Gaussian and can be simply done. Now we turn to 

the construction of ^f, in- 
One can choose local or nonlocal forms for ^i n j, the latter are obtained by insertion of 

<f)(x, y) (27) in the local expressions, which are 
ty inJ = ^Tiip for mesons, I\ = 1, 75, 7^, 7^75, 

^inj — il>Fif(a, Da)ijj for hybrids, where / is a polynomial, the simplest form is ip^ia^ip 
*&in,f = traA 2 a and £r(A 3 a 3 ) for glueballs made of 2 and 3 gluons respectively. A 2 , A 3 are 

polynomials made of D U (B) and ensure the needed tensor structure for given quantum numbers. 

In the simplest case A 2 = A 3 = 1. 

^inj = e abc4'f 1 a' l P b f2l3' l f jC f3'y^^ ia ^ 213 ^ 3 ' 1 ^ OY baryons, where a, /j, a are color, flavour and Lorentz 
index respectively. 

It is clear that higher hybrid states for mesons and baryons are obtained by adding additional 
factors of a n in ^i n j- 

Several comments are in order. Our definition of hybrids differs from that of the flux tube 
model [35] , where hybrid state is a result of the string excitation. In the latter case the quantum 
numbers of the hybrid are fixed by the mode of the string. In our case gluon a M sitting on the 
string has its own spin and makes account for larger set of valuesQ. Also dynamics is different 
and local state ^m,/ for hybrid in the flux tube model does not exist. The same differences exist 
for glueballs. In contrast to other definitions our definition of hybrids and glueballs is unique 
and is given in the field theoretical terms. The correspondence with the lattice is also clear: to 
the gluonic excitation in there corresponds a plaquette on the parallel transporter, which 
is field strength F^ in the continuum limit, (or a shift corresponding to an extra link, which is 
A a or D u in the continuum). 

Insertion of given ^m,/ in (45) yields after integration over Da u DipDi/j the hadron Green's 
function, which can be written [36] respectively for mesons 




(47) 



G q (X,Y) = {-id-im-gB)x 



(48) 




(49) 



G gluebaU = (tr(A¥\X)G g (X,Y)AT\Y)G 9 (X,Y))) B 



(50) 




G B (X,Y) = (tr(eK^Y[G^(X,Y)eK^)} 



B 



(51) 



the author is grateful to Yu.Kalashnikova for a discussion of this point 
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The structure of hadron Green's function is clear: it is a product of quark and gluon Green's 
function averaged over the background field B^. In the next section we shall discuss the way, 
how to make this averaging explicit and express it through the Wilson loop. 

So far we have neglected interaction of propagating quarks and gluons with the perturbative 
field a M , i.e. color Coulomb exchanges. For quarks those can be easily restored and summed up 
in the exponent [37,38] as we shall show below; for propagating gluons perturbative exchanges 
do not reduce to color Coulomb interaction, and should be treated perturbatively term by term 
[39]. It was discussed in [39] that next-to- leading corrections strongly damp the main term. 



5 Relativistic path integrals and relativistic Hamiltoni- 
ans 

In this section we shall discuss relativistic path integrals in the form of Feynman-Schwinger 
representation (FSR), first introduced for hadron Green's functions in [36] and later developed 
in [37], (for earlier references to FSR see [36-38] and [42]) and relativistic Hamiltonians which 
were obtained from FSR in [4,5]. These Hamiltonians will be used in the next sections to 
calculate relativistic spectrum of mesons and baryons, and also hybrids and glueballs. 

The basic approximation which will be used in this section is the perturbative treatment 
of spin degrees of freedom, i.e. of fine and hyperfine spin interactions. As one will see, for 
most mesons and baryons (except for Nambu-Goldstone mesons n and K) this is a good 
approximation and the resulting spectra are in good agreement with experiment. 

A more general treatment, valid also for Nambu-Goldstone mesons and displaying chiral 
symmetry breaking, will be discussed in the second part of these lectures. 

We start with the meson Green's function (47) and for simplicity consider flavour nonsinglet 
meson, so the last term on the r.h.s. of (47) is absent. 

For the quark Green's function in the external field = + a M one can write the FSR 
[36,37] 

G q (x, y) = (m-D) ds(Dz) xy e- K $ a (x, y) (52) 

J 

where notations used are 

2 i 1 f S J /d- Z p. \2 



K = m s + t / dr(-^Y; D, = d^- igA^ (53) 
4 jo or 

(Dz) xy = Um N ^ ft ^ e w(E«»)-(*-»)), Ne = s (54) 

rx rs 

$> a (x,y) = P A P F exp(ig A^dz^exp^g dra^F^), (55) 

J y JO 

<V = ^(7m7^ - 7^7m)> £W = z (n) - z(n - 1). 

One can write the same FSR for the antiquark Green's function; taking into account that 
in this case one should use the charge-conjugated field 

Af) = C-*A,C=-Al FW = -F%, (56) 



and 
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the ordering of PaPf in (55) changes its direction as well as limits of integration (x, y) — > (y, x). 
As a result one obtains the meson Green's function as the path integral over closed Wilson loops 
with insertion of magnetic moment terms o^F^ of different signs for quark and antiquark, as 
it should be, 

G M (x,y) = (trrW(m - D) J ds J dse~ K ~ K ' {Dz) xy {Dz) xy T^\m - D)W F ) (57) 

Here the bar sign refers to the antiquark, and 

W F = P A PFexp(ig j> dz^A^)exp(g dTa { ^F^ v )exp(-g j°G$F^df) (58) 

In (57) enter integrations over proper times s, s and r, f, which also play the role of ordering 
parameter along the trajectory, = z^(t), z^ = z^{f). 

It is convenient to go over to the actual time t = z± of the quark (or antiquark), defining 
the new quantity /i(t), which will play very important role in what follows 

2/i(t) = ^, t = z A (r) (59) 

For each quark (or antiquark and gluon) one can rewrite the path integral (52), (57) as follows 
(see appendix 1 of the second paper in [5] for details of derivation) 

J ds(D A z) xy ... = const J Dfi(t)(D 3 z) xy ... (60) 

where (D 3 z) xy has the same form as in (54) but with all 4-vectors replaced by 3-vectors, and 
the path integral Dji(t) is supplied with the proper integration measure, which is derived from 
the free motion Lagrangian. 

In general fj,(t) can be a strongly oscillating function of t due to the Zitterbewegung. In 
what follows we shall use the steepest descent method for the evaluation of the integral over 
D/j,(t)i with the extremal //o(0 playing the role of effective or constituent quark mass. We 
shall see that in all cases, where spin terms can be considered as a small perturbation, i.e. for 
majority of mesons, /j,q is positive and rather large even for vanishing quark current masses 
m, to, and the role of Zitterbewegung is small (less than 10% from the comparison to the 
light-cone Hamiltonian eigenvalues, see [40,41] for details). 

Now the kinetic terms can be rewritten using (59) as 

K+R - jf*tw + ^ [(m? + 1|} (61) 

where T = x 4 — y 4 . In the spin-dependent factors the corresponding changes are 

rs rT dl 

Jo dT ^ = /. Mf)"-"'^ (62) 

In what follows in this section we shall systematically do perturbation expansion of the spin 
terms, since otherwise (57) cannot be estimated. Therefore as the starting approximation we 
shall use the Green's functions of mesons made of spinless quarks, which amounts to neglect in 
(57), (58) terms (to — D), (to — D) and o^F^. As a result one has 

G { $(x,y) = const J Dfi(t)Dji(t)(D 3 z) xy (D 3 z) xy e~ K - R (W). (63) 
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Our next approximation is the neglect of perturbative exchanges in (W), which yields for 
large Wilson loops, R, T > T g , 

(W) B = const exp(-aS min ) (64) 
where S m i n is the minimal area inside the given trajectories of quark and antiquark, 

S m in = f o dt dp^Jdetg, g ab = d a w M <9 fe u^ ', a, b = t, (3. (65) 

The Nambu-Goto form of S m i n cannot be quantized due to the square root and one has to 
use the auxiliary field approach [42] with functions v(/3,t) and r)((3,t) to get rid of the square 
root, as it is usual in string theories. As the result the total Euclidean action becomes [5] 

T f 1 jar 1 /™ 2 , ^ \ , M*) d2 



A = K + K + aS mm = dt d(3{-(— + —) + ^^R 2 + 

Jo Jo 2 fi(t) n{t) 2 

+ + ^2 + £y r 2 _ 2r] ^ r) + ^2]} (66) 

2 2 v 

Here /i + = /x + p, p — -^r, i?j = M ^^ Zi , 1 = Zt — Zj. Note, that integrations over /i, v 
and 7] effectively amount to the replacement by their extremum values [5]. 

Performing Gaussian integrations over and r] one arrives in the standard way at the 
Hamiltonian (we take m = m for simplicity) 

TT p 2 r +m 2 L 2 /r 2 



Mr) PW p + 2!i{p-\)*v{p)dp 

a 2 r 2 ri d/J ,i 
+ —JoW) + J»— d ^ (67) 

where p 2 = (pr) 2 /r 2 , and L is the angular momentum, L = (r x p). 

The physical meaning of the terms /i(t) and can be understood when one finds their 
extremal values. E.g. when a — and L = 0, one finds from (67) 

H = 2yf^ + m 2 , /i = y/p 2 + m 2 (68) 

so that no is the energy of the quark. Similarly in the limiting case L — > oo the extremum over 
z/(/?) yields: 

"o(/?) = y = f= , H 2 =2rra^L(L + l) (69) 

so that vo is the energy density along the string with the (5 playing the role of the coordinate 
along the string. 

Let us start with the L = case. Taking extremum in v((3) one has 

Hi = + Hit) + or. (70) 

Here /i(t) is to be found also from the extremum amd is therefore an operator in Hamiltonian 
formalism. Taking extremum in /j,(t) one obtains 



H 2 = 20o 2 + m 2 + or. (71) 
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The Hamiltonian (71) is what one traditionally exploits in the relativistic quark model (RQM) 
[6-8] (apart from color Coulomb and spin-dependent terms to be discussed below). The RQM 
was an essential step in our understanding of hadronic spectra, see especially [7] for an en- 
cyclopedic survey of model predictions. At the same time the usual input of RQM contains 
too many parameters and the model was introduced rather ad hoc. Another deficiency of this 
model at this point is twofold: 

i) one usually takes m in (71) to be constituent quark mass of the order of 100-200 MeV, 
which is introduced as an input. Instead we have in (71) the current quark mass renormalized 
at the reasonable scale of 1 GeV. Hence it is almost zero for light quarks; 

ii) the form (71) is used in RQM for any L, and to this end one writes in (71) pi — > p 2 . 
However the Regge slope for both (70) and (71) is l/8a instead of the string slope 

which occurs for the total Hamiltonian (69), since the RQM Hamiltonian (71) does not take 
into account string rotation. 

Still for L = the form (71) is a good starting approximation, and it is rewarding that our 
systematic approach makes here contact with RQM. Sometimes it is convenient to use instead 
of (71)a more tractable form (70) where /i(t) = /io and /x is to be found from the extremum of 
the eigenvalue of Hamiltonian Hi. 

H^ = Mii>, Mi = Mi(/2o), ^ = 0. (72) 

C/i 

This procedure was suggested in [4], and its accuracy was tested recently in [43]. In Table 1 
listed are eigenvalues of (72) and (71) (M^ and respectively) for several values of the 
radial quantum number n. 

One can see that the difference between and M^f 1 is around 5%. 

Equation (70), (72) can be reduced to the Airy equation with the eigenvalues written as 
(for m — 0) 

M t l 1 )=4 / i (n), /io(n) = v^(^) 3/4 (73) 

where a(n) is the corresponding zero of the Airy function. The Table 2 gives a set of lowest 
a(n) and no(n) for L = 0, 1,2. The Ho(n) for nonzero L values, L — 1,2 are used later for 
spin splittings of obtained levels. The first comprehensive study of the equation (70) with fixed 
fi(t) = m was done in [44], where also a{n) are quoted (called there ((n,L)). 

Finally we turn to the Hamiltonian (67), where the eigenvalues have been found in [43] 
using WKB method. The corresponding values are listed in Table 3 for the case m — 0. 

An equivalent numerical method was used to find eigenvalues of the rotating string with 
quarks at the ends in [45], which should be compared with the entries of Table 3 and agree with 
accuracy of about 5%. Comparing mass eigenvalues in Table 1-3 one can notice that they are 
too heavy to reproduce the experimental values, and one needs to subtract some 700-800 MeV 
to get to the region of p-meson. In RQM one introduces by hand a large negative constant in 
the Hamiltonian, subtracting just this amount. In our approach here when spins are treated 
as a perturbation we do not understand how this large negative constant occurs. Indeed, there 
appears a constant Co due to the finite string width, but this is too small Co ~ 0.1 -=-0.2 GeV to 
explain the effect. It is clear however, that this constant is due to nonperturbative selfenergies 
of the quark and antiquark, and should be sensitive to the quark current mass (i.e. flavour) 
and chiral symmetry breaking. As we shall show in the second part of these lectures (to be 
published separately), the solution of this problem, which may be called" the Regge intercept 
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problem", lies in the correct treatment of spin and chiral degrees of freedom, and to obtain a 
correct Regge intercept of meson trajectory one should go beyond the perturbation theory for 
spin effects. 

Generalization of (67) to the case of unequal current masses of quarks is straightforward 

[60] 

m 2 m 2 a + p 2 L 2 /r 2 



2/i 2/i 2 2/2 2(^(1 - C) 2 + /i 2 C 2 + Jo W-OM/5)) 
<7 2 r 2 yi _d£_ i v{(5) 



where we have defined ( 



H+ J I3vdf3 
H+p,+ J vdfi ' 



For L = one can as before find extremum in v((3) and obtain Hamiltonian 

# 1 = 1 h — + — + or (75) 

2/i 2/2 2 2/i v ; 

A simplified procedure of solving (75) as in the case of equal mass, (72), reduces to the solution 
of Schroedinger-like equation and to the finding /i, p, from the extremum of the total mass [46] , 

= M(n,p,)V, M = M(fi Q ,p, ) , 

dM dM 
dfx d/io 

here M(/i, Jx) can be expressed through the same standard numbers a(n) as given in (73) and 
Table 2, namely 

MOjlP) = ^ + ^ + ^ + = (2/i)- 1 /V 2 / 3 a(n). (77) 



Extremum conditions for \x and \x look like 

m 2 1 <T 2//3 a(n) / u~ 4//3 /i 2 / 3 

~V + 2~ 2 1 /33( At + £)2/3 



= (7£ 



m 2 1 a 2 / 3 a(n)/i- 4 / 3 /i 2 / 3 _ 

2/i 2 + 2 2 1 /3 3 ( At + /i )2/3 ^ 

/,From (78), (79) one can see that whenever the current mass is large, m <C y/a, the 
constituent mass is close to the current mass, indeed first terms of expansion are 

Till now perturbative gluon exchanges and spin dependent terms in Green's functions and 
Hamiltonians have been neglected. Now we are going to restore their contribution and consider 
to this end the general form (57). 
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Leaving out for the moment spin terms a^F^, we have to do with the Wilson loop con- 
taining both perturbative and NP fields, W(B + a). For the purely perturbative interaction, 
one can use theorem of cluster expansion to get 

1 2 

— (trW(a)) = exp[- 9 - j j dz^dz' u{{a,(z)a u (z'))) + 0(g 4 )] (81) 
Introducing perturbative Green's function (in the Feynman gauge for simplicity) 

((a,(z)a u (z>))) = ^(z - z')- 2 (82) 

and integrating out z' 4 — z A in dz±dz' A = dz±d{z! A — 24), one obtains the color Coulomb interaction 
in the exponent 

^r(trW(a)) = exp[C 2 a s [ + 1 (a s ) + 2 (a 2 s )] (83) 

1 v c J Z Z I 

where 0\(a s ) is the radiative correction coming from the spacial integral dzdz' in (81), and 
02(0;^) is the corrections coming from quartic and higher correlators of the form 

((a(l)a(2)a(3)...a(n)))dz(l)...dz(n),n> 4. 

Coming now to the original Wilson loop, W(B + a), let us write it in the form 

±tr(W(B + a)) B , a = ^tr([W(B)W(a)W mt (B,a)}) B , a (84) 

where Wi n t(B,a) is the interference term, involving both fields B^ and a M . Part of this inter- 
ference was taken into account in the previous chapter, when the freezing of a s was considered. 
This shows that interference may change behaviour of perturbative terms at large distances. 
Another example of interference was considered recently in [30,31] when treating the small dis- 
tance behaviour of static potential. In both cases it is seen that the main domain of perturbative 
fields is small distance domain, while that of NP fields B^ is the large distance domain. From 
the static QQ potential one can deduce that both regions join smoothly at r c pa 0.25/m and 
perturbative force is strong for r < r c (and stronger than both NP and interference), while NP 
force is strong for r > r c and much stronger than both Coulomb force and interference. So our 
simplifying assumption is that one can disregard interference term Wi nt altogether, and take it 
later as a correction. 

As a result confining term from W(B) and color Coulomb term enter as a sum, hence it is 
legitimate to take color Coulomb term as an additional term to be added to the Hamiltonian 
(74) or (75). Note that this is true both for heavy quarks (nonrelativistic situation) and light 
quarks (relativistic situation). What is still left over, apart from interference, is radiative 
corrections 0(a s ) and all higher terms in 0{a 2 s ). The reason why the Coulomb term is retained 
in (83) as compared to those suppressed, is that Coulomb term is enhanced as 0(1 /v) for slow 
quarks and the quadratic term in (81) generates Sudakov-type asymptotics [38]. 

As was mentioned above, color Coulomb interaction for the qq system (situaton will be 
different for gluon systems, see below) is exponentiated in (83) which implies that Coulomb 
term can be taken to all orders, i.e. one should consider the Hamiltonian 

H c = H' — f^l (85) 
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where H' is the string Hamiltonian (74) or (75). 

One can express the eigenvalue e n (p) of the H c through the eigenvalue of the reduced 
equation, as was done before without Coulomb in (73). 

e n (/2) = (2/i)- 1 /V/ 3 a(A,L,n) (86) 

and the reduced equation looks like [45, 46] 

{ ~Jp~ 2 + |p| " R Mp) = a(A ' L ' n)xip) (87) 

where 

4a s (2/}) 2/3 p 
X = 3* o-V3 ' r = (2/2) V3 ' a( - A = °' L ' n ^ = °^ L ' ^ 

The eigenvalues a(X,L,n) are given in Table 4 for L = n = 0. 

Finding as before /x and fi from the extremum condition on M n , 

m 2 a rh 2 R 
M n = — + S + ^ + ? + £ "(/i), (89) 



one obtains equations [46] 



2/i 2 2/i 2 2 



^ + ^^(a(A) + fe |2A) (90) 



2 



(2/i) 2 / 3 <7 2 / 3 , ,9a, 



- ("W + laA |2A) (91) 

one should note, that Coulomb attraction strongly affects the values of p, p and wave function 
at origin. 

Connection between eigenfunction (p n (r) of the Hamiltonian (85) and eigenfunction Xn(p) 
of the reduced equation (87) is 

V? n (r) = (2fla)^Xn(p)Y Lm {e,ip). (92) 

For purely linear potential, i.e. for A = 0, the value of the function (p n (0) does not depend on 
n [45,47] 

l^ 0) (0)| 2 = ^, XST ) (0) = 1. 03) 

Taking into Coulomb potential, i.e. A 7^ 0; one gets an enhanced value of Xn(0), which is given 
in Table 4. Using these entries and (92) one can calculate |y? n (0)| 2 for any qq system, which will 
be used in the next chapter for the hyperfme interaction and for the quark decay constants. 

For this quantity the asymptotic freedom is also important and should be taken into account. 

The coefficient p(AF) = | x x ff | 2 , where x as (p) is the solution of (87) with ct s (r) in A de- 
creasing according to asymptotic freedom from the maximal value with Aq C d = 140MeV 
is given in Table 5. 

Now we turn to the spin terms, considering them as perturbation (a complete treatment of 
spin terms within new formalism will be given in the last section) In the derivation of spin- 
dependent interaction valid also for massless quarks the advantage of the present method will 
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be evident, since the usual treatment exploits expansion in inverse powers of quark masses, 
appropriate for heavy quarkonia only. In our case the problem reduces to the calculation of 
averages in (57), (58), with the use of (62). This can be written as [48,49] 

w - mf"v^m )ex '' { - I ik^^m m) (94) 

where for (W) one can use the cluster expansion (28). 
In the Gaussian approximation one has (Pp — > P a ) 

(W F ) = P^expiK, + K 2 + K n + K 22 + K 12 ) (W) (95) 

and e.g. 

K> = tg 2 1 ^y<T« / d S ^(F,Mt))F>. p («')), (96) 

K » = ^ W) [ (W^Mm, (97) 

while other terms are obtained by replacing a^J by affl o r vice versa. In what follows 
we do not need the selfenergy terms K n ,K 22 ~ (o"^) 2 - 

Now one use representation (33) to express field correlators through D(u), D\(u). To define 
the minimal surface for integration in (W), (Wf) one connects trajectories of quark z(t) and 
antiquark z(t) by a smooth line w(t,/3) at t — i, < (5 < 1. One can approximate it by a 
straight line for simplicity: w(t,(3) = z(t)/3 + z(t)(l — (3). 

Introducing angular momenta in Minkowski space 

= (r x pj)i = ie ikm r k (t)/iz m (t), (98) 
Lf ] = (r x p 2 )i = ie ikm r k (t)pz m (t), (99) 

one has 

1 (-P)L^ 

ds ik = {dwiW k - dw k Wi) = -d(3dte ikm ( — — H ^- nL -) (100) 

i /i /i 

We refer the reader for details to [48-50]. 

The resulting spin terms can be written in the Eichten-Feinberg notations [51], namely 

T/ . . M^lV af ) Lf\.lde 2dV 1 , 



2/1/1 r dr 12/i/x 



V,{r)+ (101) 



' -^n^n k -^af)V,{r) 



12/i/i 

and Vi(r),i = 0, 1, 2, 3, 4 are expressed through D, D\ [48,50] as follows. 



71F = -L iv L t^-tW^ (102) 
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V A = / dv(3D(r,v) +3D 1 (r,iy) + 2r 2 — ^) (105) 

J —oo OT 

lde(r) f°° , [ r d\ T ^.^ . „ . ,, 9 o,dD x , 

r^r = L dV L V lOM + D lM + ^ + ^] (106) 

Equations (93-97) give the NP spin- dependent interaction through functions D,Di; the latter 
have been measured on the lattice [22 - 24] to be of the form 

D(x) = D(0)exp(-^), D x {x) = J D 1 (0)e^(-M) (107) 

with §^ = and T 9 = 0.2(0.33)/m for quenched [22,23] (dynamical [24]) quarks; while 

D(0) = 0.073Ge\/ 4 for quenched case [22]. 

Insertion of (107) into (102-106) yields an estimate of all 5 potentials e(r), Vi(r),i = 1, 2, ...4. 
The first, e(r) is just scalar potential V(r) appearing in (37), and present in (75) in the 
form or; Vs(r) and V^r) are numerically small as compared to the perturbative contribu- 
tions to be discussed below. The largest contribution comes from the spin-orbit term; V so ; 
since \^ 1 / (r)| r _ >00 = —^'( r )\r^oo = —< J , one obtains for the equal mass case the asymptotics 

(if = -if = U) 

V„(r) = "^/M, /(r-oo) = l. (108) 

For small r the behaviour of /(r) is changed, /(0) = 0, and Ko( r — > o) = const. For more 
details and discussion of the behaviour of i>j(r) and comparison with lattice data see [50]. 

Analysis of spectra of heavy quarkonia made in [32] strongly prefer /(0) = 1, in contradiction 
with (101-103). The way out may be seen in the interference terms at small r, studied in [30,31], 
which show that linear dependence of e(r) in r may persist at very small r in agreement with 
a recent dedicated lattice analysis [14]. 

We turn now to the perturbative contribution to spin-dependent interaction and to this end 
one can use again the general expressions (101-106) where instead of D,Di one should insert 
the corresponding perturbative expressions. To the order 0(g 2 ) one has 

1 6ry 

£>(*) = 0, £>!(*) = g-J (109) 

Insertion of (109) into (101-106) yields 

V 1P = 0; ±V> p (r) = g|, V 3p = ^, e p (R) = V, P (R) = ^p^(R). (110) 

The final expression for the mass of the meson with spin interaction taken into account pertur- 
batively is obtained from (75) and (85), (101) and has the form 

M ra = M ra (^ , / , ) + 32SlS2TO > 2 (0) + 
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2/X5 2 ^ r /i /i 2/xg 3 r 3 



(3sins 2 n - SiS 2 ) 4 a 5 



Where S = Si + s 2 , and 



Mn(/io, ^o) = ^ + v + ^ + V + + AM - C = + AM - Co (112) 

Here AM takes into account the so-called string correction, i.e. the difference of eigenvalues 
of the Hamiltonian (67) and (70) due to the string rotation. According to [5] this can be written 

as 

AM = _ W + i) . (113) 

3(M,< 0) )3 

Let us note the main properties of the mass formula (111). First of all it has the minimal possible 
number of input parameters; those are current quark masses (defined at the appropriate scale 
of lGeV), string tension, a s (r) or equivalently Aq C d and the overall subtraction constant C . 

We do not introduce constituent quark masses and do not expand in inverse quark 
masses to get spin interaction, as it is common in relativistic quark models. In addition we 
take into account string rotation, which is not done in those models. 

Moreover, our expression of the total mass M n is the result of derivation from the first 
principles of QCD with all steps and approximations clearly visible, so that one can check and 
improve if need be. This is in contrast to the ad hoc model building principle. 

Now we come to the confronting of our approach and first of all mass formula (111) with 
experiment and lattice data. 



6 Spectrum of light mesons, hybrids and glueballs 
6.1 

The spectrum of mesons and glueballs (and of any two constituents connected by the string 
with spins and Coulomb force neglected) is given by the Hamiltonian (74). 

The denominator in (74) is actually the moment of inertia of the system, comprising both 
quark (the p term) and string (the v term) contributions. 

The eigenvalues of (74) have been found quasiclassically in [44] and given in Table 3 for the 
case of zero-mass quarks. For the case of two-gluon glueballs the corresponding eigenvalues are 
obtained by multiplication with the square root of Casimir factors ratio; since for the adjoint 
string a adj = & fund c 2 f fund) = \ a fund In Fig. 1 shown are trajectories of eigenvalues of (74) as 
functions of L and n. The Regge slope is close (within 3%) to the string slope {2i:a)~ 1 . Note 
also the linear trajectories in n with approximately 2 times smaller slope. 

The resulting spectrum for light-light mesons was found first in [34] and is close in basic 
features to that of Isgur and Godfrey [7] for L = 0, but differs for L > 0, since we take 
into account string rotation, and differs somewhat in spin splittings, since we have used for 
spin-dependent (SD) terms expressions derived from QCD (102-106), (110). In Fig.l also 
shown is the experimental p type Regge trajectory, which has almost the same slope, but 
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shifted downwards. For this reason we introduce at this point (as in most models) the negative 
constant C , as a surrogate of the nonperturbative selfenergy term for each quark, while we do 
not introduce such term for the valence gluon. The reason for this procedure and an estimate of 
this selfenergy will be given in the second part of lectures (published separately), where effects 
of chiral symmetry breaking are taken into account. 

An example of the calculated spectrum of light-light mesons is given in Table 7 for I — 1 
case and compared to experiment. It is seen that the overall agreement is reasonable and 
comparable to that of Isgur-Godfrey model [7]. Note, however, that in [7] there are 14 fitting 
parameters, while we have in our approach only 3: a s , a and Cq. The most important fact is 
that spin-dependent potentials are exactly derived with masses fio,fto without any expansion 
in inverse mass powers assumed in all potential quark models; moreover /io,/^o are computed 
and expressed through \fo. 

Several words should be said about pion. The present approach as a whole is based on the 
assumption that spin corrections can be computed as a perturbation. This also implies smallness 
of chiral symmetry breaking effects. Now for the pion this is not valid, hyperfine correction 
is too large, and what is even worse, if one tries to treat the spin-spin term — ^f5^(f) to 
higher orders, then the minimization of the mass M n (111) in /i is no more possible, since at 
fi — > the spin-spin term diverges. This happens only in the pionic channel and signals the 
instability of the vacuum and the onset of chiral symmetry breaking. Therefore pion should 
not be considered in the framework of the approach described. 

For radially excited pion states the spin-spin correction is smaller and the present approach 
may be more reasonable. Some of these states are presented in Table 7. 

6.2 Hybrids 

Hybrids and glueballs are the most interesting examples where all advantages of our approach 
can be seen clearly. First of all this is the definite construction of the hybrid ( and glueball) 
Green's function (49) based on the background perturbation theory. Note the difference of our 
definition of hybrids from that of the flux-tube model; at the same time the structure of in 
and out hybrid states ^m,/ is similar to that of the lattice QCD. Moreover, since for hybrids 
we do not introduce any new parameters (and effective mass of the gluon ji g will be calculated 
through string tension), one can predict hybrid and glueball masses unambiguously, in this way 
seriously checking the whole formalism. The first studies of hybrid spectra in this method have 
been done in [52-54]. For later development using variational methods see [55]. 

Let us turn to the construction of the hamiltonian for the hybrid system. We start with 
the hybrid Green's function and treating spin terms as perturbation, disregard in the first 
approximation the gluon spin term 2igF fMU in (9). Using FSR we get similarly to (57) the form 

^ /*oo /*oo /*oo 

G hyb (x,y) = (trT«\m-D)j ds J ds J ds g e- K - K - K ^{Dz) xy {Dz) xy {Du) xy 

xT^(m-b)W hyb ) (114) 

where 

W hyb = ($(x, y)) aP t a ^ ab (x, y)t b Sa <5> lS (y, x), (115) 

and y), (<&(?/, x)) are transporters (27) belonging to the quark (antiquark) Green's function, 
while $(x, y) is that of the gluon in the adjoint representation. At this point one can use the 
large N c approximation, which yields accuracy around 10% in the most problems of QCD. 
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Then the gluon line can be replaced by the 'tHooft rule with a double qq line, 

tfr&abtsa -> Qfiaix, V^iV, x) (116) 

and the vacuum average of Wh y b reduces to the product 

(W hyb ) = (W(C 1 ))(W(C 2 )) (117) 

where C\, C 2 are adjacent contours made of quark (antiquark) trajectory and gluon trajectory. 

The subsequent transformations with ( |114j ) are the same as in the qq case leading from (57) 
to (67). 

We shall use only the simplest form of the Hamiltonian for zero angular momenta, equivalent 
to (70) and (112). In the case of the hybrid state it is 

m\ m\ Mi + /i 2 + us . , n ( - \ / 11s s 

H hyb = — + — + +h H - C {qqg) (118) 

with 

h «-^w + w )+ ^ m -^ l+ ^ w -^ ; (119) 

where we have defined c.o.m. coordinate and total effective mass respectively as 

i* = £^-, /i = I> (120) 

k=l A* k=l 

while Jacobi coordinates £,i,r]i, are expressed through those of q, q: zf^ and of gluon zf^ as 



V A*A*3 



^) = ^ + / m ^ 6 + / m ^ 2 

y fjL(fJLi + fi 2 ) % V A*i(a*i + /"a) 



z (2) = ^ + / m /^3 ^ _ / ^1 , 121 x 

y/i(/ii + /i 2 ) 4 Y /x 2 (//i + // 2 ) 

Here m is an arbitrary mass parameter. At this point one should comment on gluon ex- 
changes between the valence gluon in the hybrid and quark (or antiquark). As we shall see in 
the next section, devoted to glueballs, a valence gluon does not have a color Coulomb interac- 
tion with another gluon or quark: the OGE diagrams existing in the lowest order do not sum 
up in the ladder-type series of the Coulomb interaction, since exchanged gluon and valence 
gluon are identical and moreover four-gluon vertex is also important. This is in contrast to 
the qq or qqq systems, where the Coulomb series occurs naturally as a first term in the cluster 
expansion series, see (81). The perturbative series for the gluon-gluon case is summed up in 
the Lipatov approach [56] and the result for the gluon ladder is strongly damped by the one- 
loop corrections. Therefore we disregard here the corresponding Coulomb terms for qg and qg 
systems. 

The Hamiltonian (118) can be used in problems of two types. First, one can consider, as it 
is done on the lattice, [57], fixed (very heavy) quark and antiquark at a distance R from each 
other, and gluon in some angular momentum and spin state. 
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In this case one can define as in [57] the potentials V r (r) for the static QQ pair at distance 
r with gluon in a state, which can be denoted as in the diatomic molecule, i.e. with projection 
of the total gluon momentum J on the axis r, 

A = Jr/r, A = 0, 1, 2, ...(£, II, A, $) (122) 

and combined operator of inversion w.r.t. midpoint of QQ and charge conjugation i] cp = 1(— 1) 
denoted as g(u). £ can be also even (odd) w.r.t reflection in the plane of QQ, denoted as +(— )• 
The Hamiltonian for the gluon in such two-center problem looks like 

H QQg - H QQ = f + |1 + alzW - z Q | + _ z ^| + y L5 + 1?± (123) 

where the spin-orbit term is similar to that in (101), while the Coulomb term is a result of 
subtraction from the coulor-octet repulsion ^ in the QQg system of the attractive color-singlet 
term —4^. 

The Hamiltonian (123) can be studied numerically. 

A rough estimate of the eigenvalues of (123) is obtained when one expands string potentials 
in (123) assuming that gluon slightly vibrates around the middlepoint of QQ, expanding the 
square roots. 

Then the spectrum is 

AW^) = (v) 1/3 (-+^) 2/3 + ^^ n = 0,l,2... (124) 

A comparison of (124) with the lattice calculation [57] shows a good qualitative agreement. 

We turn now to the real physical objects, qqg systems with light or heavy quarks. In this 
case one should take quark kinetic energies into account and consider qqg as a 3-body system 
with the Hamiltonian given in (118). The most efficient way to treat this problem is the 
hyperspherical expansion [58]. One defines the hyperradius p as 

p 2 = V 2 + e, (125) 
and expands the wave function in a series 

*(tl) = ExM(^ (126) 

Kvn 

where Qku is hyperspherical angular function [58], and v denotes all quantum numbers in 
addition to the global momentum K = 0, 1, 2, ... 

Keeping only one harmonics with the given K, one obtains radial equation (neglecting for 
the moment Coulomb and spin-dependent terms) 

d 2 

["2^2 + W (P)]XKM = SKnXKnip) (127) 

where 

W(p) = + ap—(a 13 + a 23 ) (128) 

2mp z 157r 
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and 

a = if+-, a i3 = J , 2 = 1,2. (129) 

2 V 

For an estimate with some 5% accuracy, one can find Ek n by simply minimizing W(p) in p. 

Assuming now equal quark masses, mi = m 2 one could obtain pi = p 2 , and taking the total 
hybrid mass as 

2 

M x „ = — + p x + ^ + e Kn - C (qqg) (130) 
Pi i 

and finding the minimal value Sku from W(p ), one has 

M Kn = ^ + /ix + ^ + - C te) (131) 

where (7 = (^) 2 / 3 [q;(q; + l)] 1 / 3 . It is reasonable to take C (qqg) the same as in the qq case 
since the constant term is due to quark selfenergy terms only. 

In the limit of heavy mass mi ^> y/a, one get from minimization pi ps m 1; and minimizing 
in p 3 one obtains 

M^ n = 2mi + 2v^C 3 / 4 - Co(g^) (132) 
In particular, subtracting the heavy mass 2m! — C one gets 

AM 9/^/4 2.72V? = 1.375GeV, ^ = 

For light quarks one takes the same C for a light hybrid as for the corresponding light meson. 
For the lowest exotic hybrid 1 h one should take K — 1 (L = 1 for the gluon) and agreement 
with recent lattice calculation of the bbg mass excitation is very good, see Table 8. 

6.3 Glueballs 

This section is based on papers [59,60]. The L = Hamiltonian for glueballs (neglecting spin 
and perturbative interaction) is obtained from that of qq system by replacing <Jf un d by a a dj- 

H' Q = — + p + a adj r (134) 

The value of o a dj in (134) can be found from the string tension of qq system, since the Casimir 
scaling found on the lattice [33] predicts that 

C 2 (arfj) 9 
aadj ~ C 2 (fundf fund ~ 4 a/W (135) 

For light quarks the value of crf un d is found from the slope of meson Regge trajectories and is 
equal to 

<Jfund = ^7 ~ O.lSGeV 2 (136) 

/,From that we find 

a adj ps OAOGeV 2 (137) 
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In what follows the parameter \x and its optimal value /j,q, which enters in (134) play very 
important role. The way they enter spin corrections in (101) and magnetic moments shows 
that fiQ plays the role of effective (constituent) gluon mass (or constituent quark mass in the 
equation for the qq system) . 

In contrast to the potential models, where the constituent mass of gluons and quarks is 
introduced as the fixed input parameter in addition to the parameters of the potential, in our 
approach /x is calculated from the extremum of the eigenvalue of equation (134), which yields 

A*o(«) = v^(^) 3/4 ,M (n) = 4// (n) 

where a = a a dj for gluons and a = 07 for masses quarks, and a(n) is the eigenvalue of the 
reduced equation ip + (a(n) — p — L(L + 1)/ p 2 )ip = The first several values of a(n) and po(n) 
are given in the Table 2, and will be used below. 

Note that our lowest "constituent gluon mass" po(n = L = 0) = 0.528GeV (for Of = 
O.lSGeV^) is rather close to the values introduced in the potential models, the drastic difference 
is that hq depends on n, L and grows for higher states. 

/From Table 1,2,3 and Fig. 1 one can see that mass spectra of the Hamiltonian (67) are 
described with a good accuracy by a very simple formula (a similar conclusion follows from 
calculations in [45]) 

M 2 

- — — L + 2n r + c\ (138) 

27TCT 

where L is the orbital momentum, n r -radial quantum number and C\ is a constant w 1.5. It 
describes an infinite set of linear Regge-trajectories shifted by 2n r from the leading one (n r = 0). 
The only difference between light quarks and gluons is the value of a, which determines the 
mass scale. 

Thus the lowest glueball state with L = 0,n r = according to eqs. (6), (5) has M 2 = 
AMGeV 2 . 

It corresponds to a degenerate ++ and 2 ++ state. 

M = 2.01GeV (139) 

In order to compare our results with the corresponding lattice calculations [61-63] it is 
convenient to consider the quantity M/^/aJ, which is not sensitive to the choice of string 
tension aQ. The spectrum of glueball states obtained in lattice calculations is given in the 
Table 9, where masses of glueballs for values of 07 used in these calculations are given. 

^From these data we have for L = 0, n r = states the spin averaged mass 

M M(0++) + 2M(2+) 1 



a f 3 V°7 



the value 4.61 ±0.1, which should be compared to our prediction M theor (L = 0,n r = 0)/ y/ /oJ = 
4.60. 

For radially excited states our theory predicts 

jyjtheor 

■{L = 0, n = 1) = 7.0 (141) 



2 Note that the value 07 ~ 0.2lGeV 2 used in lattice calculations differs by about 20% from the "experimental" 
value <j } = O.lSGeV 2 . 
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Lattice data [61] give for thise quantity 



M lat 

: (L = 0,n r = 1) = 6.8 ±0.35 (142) 



For L — 1, S — 1 states one can define spin- averaged mass in a similar way 

M M(0-+) + 2M(2~+) 1 



lattice data [61-63] yield 



(143) 



M lat M lat 

: (L = l,n = 0) = 6.36 ±0.5; (L — l,n= 1) = 8.1 ± 0.5 (144) 



which is in a reasonable agreement with our prediction 

j^theor j^theor 

.(L = 1, n r = 0) = 5.95; (L = 1, n r = 1) = 8.0 (145) 



/cr f V a f 
For L = 2, n r = the spin averaged state has 

M theor (L = 2,n r = 0) , . 
i — = 7.0 146 

Lattice data [61] exist only for 3 ++ state which yields ~|=~^ — 7.7. The overall comparison 
of spin- averaged masses computed by us and on the lattice shows a striking agreement. 

Thus we come to the conclusion that the spin-averaged masses obtained from purely con- 
fining force with relativistic kinematics for valence gluons are in a good correspondence with 
lattice data, which implies that the dominant part of glueball dynamics is due to the QCD 
strings. 



6.4 Spin splittings of glueball masses 

Here we shall treat spin effects in a perturbative way, in the same manner, as it is done with 
spin effects in heavy and light quarkonia; a glance at the lattice results given in Table 9 tells 
that spin splittings in glueball states apart from 2 ++ — ++ amount to less than 10-15% of the 
total mass, and hence perturbative treatment is justified to this accuracy level. 
The two-gluon mass operator can be written as 

M = M (n, L) + SLM SL + S (1) S (2) M 55 + M T , (147) 

where M is the eigenvalue of the Hamiltonian H' = H Q + AH pert , and H is given in (134), 
while AH pert is due to perturbative gluon exchanges and as discussed in [60], gives a small 
correction, which will be omitted here. 

To obtain three other terms in (147) one should consider averaging of the operators F in 
(9) which enter in exponent for the Green's function and take into account that 

- 2iF uv = 2(S«B« + S (1) E«)^ (148) 
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and similarly for the term in the integral / Fdr' , with the replacement of indices 1 — ► 2. Here 
gluon spin operators are introduced, e.g. 

(^rn)ik = ~^ e miki h k = 1,2,3, {S^)i4 = —%bi m (149) 

Background gauge condition (46) allows to exclude // = v — 4, and hence S^ 1 ) in (148). Since 
the structure of the term F is the same as in case of heavy quark with the replacement of the 
heavy quark mass by the effective gluon parameter /io(see (134)), one can use the spin analysis 
of heavy quarkonia done in [46,50], to represent the spin-dependent part of the Hamiltonian in 
the form similar to that of Eichten and Feinberg [51] 

SL,2dVi 2dV 2 , S (1) S (2) T ,, . 



+ -^(3(S (1) n)(S (2) n) - S (1) S (2) )V 3 (r) + AV (150) 

where S = + S*- 2 **, AV^ contains higher cumulant contributions (which can be estimated to 
be of the order of 10% of the main term in (150)). 

The functions Vi(r) are the same as for heavy quarkonia [48,50] except that Casimir oper- 
ator is that of adjoint charges, the corresponding expressions of Vi{r) in terms of correlators 
D(x), Di(x), are given in (46,50). Both D and D\ are measured on the lattice [22-25] and Di 
is found to be much smaller than D. Therefore one can neglect the nonperturbative part of 
Vz(r), while that of V4 turns out to be also small numerically, Mss(nonpert.) < 30MeV^, and 
we shall also neglect it. 

The only sizable spin-dependent nonperturbative contribution comes from the term ^ 
(Thomas precession) and can be written asymptotically at large r as 

AH(Thomas) = (151) 
1 ' r 2/i 2 { } 

Now we come to the point of perturbative contributions to spin splittings. 

In heavy and light quarkonia these contributions are simply obtained from (102-106) sub- 
stituting for D, Di their perturbative expressions to the order 0(a s ) (109). In case of glueballs 
one should make replacement C2(fund) — > C^ifldj) in (110) and in addition take into account 
two corrections due to the fact that i) valence and exchanged gluons are identical ii) there is 
4-gluon vertex. These corrections have been taken into account in [64] and amount to some 
reduction of coefficients in (110). 



The resulting matrix elements in ( 147 ) look like 



, Avert) _ 3C 2 (adj) a s 

MsL ~ 4/, 2 V } (152) 



ii4r t} = [ r m 3) w) (153) 



(pert) h-KC 2 {adj) 
■ss 

(pert) _ C 2 (adj) a 



3/i 



M ivert) = m^ ( ^ (3S ( 1)nS(2)n _ S (1) S ( 2))) (154) 

Ho r A 
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^From ( |153j ) one can see that Mss can be written as 

M SS = -i\R(0)\ 2 (155) 

To make simple estimates, we shall neglect first the interaction due to perturbative gluon 
exchanges between valent gluons. Indeed it was shown in [39] that this interaction cannot 
be written as Coulomb potential between adjoint charges, and comparison to perturbative 
Lipatov pomeron theory [56] shows that it is much weaker than Coulomb potential. Neglecting 
this interaction altogether, one gets the lower bound of spin-dependent effects, since all matrix 
elements, like (S^ 3 '(r)), (-), (\) are enhanced by attractive Coulomb interaction. 

For purely linear potential one has simple relation, not depending on radial quantum number 
n r [44,47] 

|*(0)l 2 = |jR( ° )|2 = ^ V '^ = ^L ( i5 6) 

47T 47T 47T 



Using (72,73), one obtains 

M.. = * 

Mr 



M ss = (157) 



and for n r = 0, 1 and a s = 0.3 one obtains 

M ss (n r = 0) = 0.3GeV, M ss (n r . = 1) = 0.20GeV (158) 

For M(0 ++ ) and M(2 ++ ) one has the values given for the sake of comparison with lattice 
calculations in Table 9 for Of = 0.23GeV^ 2 and a s = 0.2(0.3) 

For L > one needs to compute spin corrections Msl and Mr. First of all one can simlify 
matter using the equation (it is derived in the same way, as (156) was derived in [47]) 

L(L + l)(l> = ^<V(r)> (159) 

For V(r) = o a ^r both M^l and Mj? ert ^ are easily calculated and used in Table 9. 

The estimate of AH(Thomas) = 1 r 2 ) is more cumbersome since this matrix element 

is very sensitive to the behaviour of Vi(r) at small r, where the asymptotics (151) is not yet 
achieved. Therefore one has to use explicit expressions (102, 103) of V[, V% through correlators 
D,D 1 (see [50] for details and discussion). 

The resulting figures for AM t homas are used in Table 9. Combining all corrections and values 
of Mo from Table 2 one obtains glueball masses shown and compared with lattice data in Table 
9 for e f = 0.23GeV 2 . 

The general feature of spin-dependent contribution AH S is that it dies out fast with the 
growing orbital or radial number, which can be seen in the appearance of the /Xq factor in the 
denominator of ( |150| ). 

Indeed, from (73) one can derive that M w 4/i and therefore AH S ~ M\ni) (P(\))-> wriere 
O stands for terms like const. - or const' .\ (from perturbation theory). Hence spin splittings of 
the radial recurrence of states ++ ,2 ++ or h ,2 h should be smaller than the corresponding 
ground states. This feature is well supported by the lattice data in Table 9. 

We end this section with the discussion of heavy-light mesons, which can be calculated 
also in the present approach, see [46]. While in [3] is given another and more exact approach 
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for heavy- light mesons, taking into account chiral symmetry breaking, here we shall only 
demonstrate that the present approach enables one not only reproduce the spectrum, but also 
calculate more delicate characteristics, e.g. decay constants f m , M = B, B s , D, D s and lepton 
width. These are characteristics computed in [46] and shown in Table 10. One can see a 
reasonable agreement of all entries with lattice data and experiment. 

We do not touch in these lectures the topic of heavy quarkonia, which has been also ex- 
tensively studied in the framework of the present approach [19,32], since this is the subject of 
lectures of F.J.Yndurain at this School [66]. 



7 Light— cone Hamiltonian and spectra for mesons. 

In this chapter we follow mostly the work done in [41]. The light-cone description of the hadron 
wave functions is widely used, since it allows to get a direct connection to the parton model and 
its QCD improvements p? |. The latter however are mostly perturbative and nonperturbative 



contributions are introduced via OPE and QCD sum-rules. In this way the concept of the string 
- the main nonperturbative QCD phenomenon - is totally lost. On physical grounds it seems 
that the string is the essential ingredient of the dynamics for large (r > 0.3 fm) distances, 
and it would be interesting to understand its contribution to the light-cone wave functions, 
formfactors, structure functions etc. 

In particular, what is the QCD string in the parton language? Should one associate it with 
the gluon contribution as an assembly of gluons compressed inside the string - or with the 
constituent quark mass? 



There two contrasting points of view have been proposed already decades ago [58, 69]. In 



|B~8H the sea quarks and gluons enter as separate entities and one could associate the gluon 



distribution with the string in the same way as photons with the Coulomb field of the charge in 
the Williams- Weizsacker method. In contrast to that in |)9| the quarks have been considered 



as constituents with structure, and string does not appear separately. Recently a quantitative 
analysis was performed [|7(J] of quark distribution in the pion starting from that in the nucleon 



and assuming the same internal structure of quarks in nucleon and in pion. It is still an open 
question how the structure of the constituent quarks is formed, and to which extent it can 
be explained by the adjacent piece of the string. This problem can be elucidated partly in 
the present approach, since our light-cone Hamiltonian contains the string on the light cone 
explicitly. It allows to separate the contribution of the string to the parton distribution, in 
particular to the momentum sum-rule. There is another source of structure in the constituent 
quark - chiral symmetry breaking which creates the chiral part of constituent quark. This 
problem will not be discussed below, see e.g. 0. 

Let us make a few comments about the reasons why the light-cone dynamics is interesting 
to explore in QCD. One simple reason was stated in the beginning: it is to formulate in terms 
of parton language and therefore to single out the new nonperturbative contribution in the 
familiar parton picture. Second reason is quite general: the necessity of using Hamiltonians 
in moving frames to calculate observables, like form-factors or cross sections, where overlap 
integrals of eigenfunctions in different frames enter. An important advantage of the light-cone 
wave functions is that they allow to calculate formfactors and structure functions directly, 
without additional boosts, which typically are the dynamical ones, i.e. require the use of the 
exact Hamiltonian. 
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There is still another and deeper reason for using the light-cone formalism. It is based on 



the expectation [71] that the field theory in general has a more simple vacuum structure on 
the light-cone and in addition some classes of diagrams (containing backward-in-time motion) 
are suppressed. Hence one may expect that on the light-cone it is possible to choose the field 
degrees of freedom which are independent and simplest. 

We consider the relativistic quark-antiquark pair with the masses mi and 7712 connected by 
the straight-line Nambu-Goto string with the string tension o in 3+1 dimensional space-time. 

Given a qq Green's function in the coordinate space G(xx;yy), where xx(yy) are final 
(initial) 4-coordinates of quark and antiquark, one can define the Hamiltonian H through the 
equation (in the euclidean space-time) 



dG 
df 



— — = -HG 



(160) 



where T is an evolution parameter corresponding to some choice of a hypersurface S. In a 
particular case of the cm. Hamiltonian the role of T is played by the center-of-mass euclidean 
time coordinate T = (x 4 + x 4 )/2 and the hypersurface J2 is a hyperplane £4 = x 4 = const. 
With the notations for the vectors a M , 6 M 

ab = CLubu = a-ibi — a &o = a±b± + a + 6„ + a_5 + , 



a± 



a-3 ± Qq 
V2 ' 



one can define the hypersurface J] through the qq coordinates z^.z^ as 



and the kinetic part of the action A 



A = K + K + aS mir . 



has the form 



K + K 



, /q i l( T "> dT + \j z l( T ) dr + L mldr+ I m 2 2 dr 











dz 4 



where we have defined 



2 v ± J 2 v ± J 2/i! 2/x 2 



2MiO+) = it- ! 2// 2 (z + ) 



or or 

The minimal surface S min is formed by connecting z^(z + ) and z^(z + ) with the same value 
of the evolution parameter z + , i.e. 

S min = a f T dz + fdPiwV - (tW) 2 ] 1/2 

JO JO 



where 



w^{z+\ 13) = z^{z + )f3 + z M (*+)(l - 13) 
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and dot and prime denote derivatives in z + and (3 respectively. 
We now introduce " center-of-mass" and relative coordinates, 



Rp = xz^ + (1 - x)'z ll , r> 



where the variable x is defined from the requirement that the term r±R± should be absent in 
the action. This yields: 

x= K + fufidp 
Mi + M2 + J vdp 

Then for the Green function one obtains: 

—A 



G(xx;yy) = J D^i 2 {z + )Dv DR^ Dr^ 



e 



where the action A 

,2 ™2 



- ( r_ + R±r ± + (((3) — x)r±r ± ) 2 



A = ^ dz+{ — + — + a^Rl + 2R_) + ar\ + 
I J I Mi A*2 



+ f-d(3rl- 
J v 



rl(Judp)- 1 
r 2 . 



(162) 



The following notation was used: 

ai = Mi + M2 + (163) 
Jo 

a = mi(1 - xf + /i 2 x 2 + f 1 v{(3){(3 - x) 2 d(3 (164) 

Jo 

Integration over DR^ leads to an important constraint: 

ai = P+ 

Furthermore we go over into the minkowskian space, which means that 

Hi -> -inf , v -> -2z/ M , A -> 
For the minkowskian action we obtain (omitting from now on the superscript M everywhere) 



To complete the Hamiltonian formulation of our problem we define canonical momenta for 
the coordinates R- and r_. As it was shown, canonically conjugated momentum to the R- is 
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P + — a±. Situation with r_ is more subtle. In order to clarify the situation let us start with 
the general form of the qq Green's function in the Feynman-Schwinger formalism 

G(x,y) = j ds dsDzDz l ~ K - R (W{C)) 

to impose boundary conditions, one can rewrite DzDz using discretization 

£ n = z(n) - z{n - 1), £ n = z(n) - z(n - 1) 
DzDz = J| d^ n dl n i dp dp'x 

n,n' 

xexp \ipQ2 Cn + X-Y) + L> + X - Y)} 

One can introduce the total and relative momenta 

P = p + P, q = 



2 ' 

and R = AR n /e; Ne = T, one has 

£,n%n ~l" £n(l x n) AR n 

Expressing £„, £ n through AR n , Ar n and going over to the momentum representation of G one 
obtains 

G{P) = fdqJl dAR n dAr n ,x 

n,n' 

xexp\iA + iP^AR n + i$^.P(l -2a;) +g)Ar„l 

where A = Jq drC = J T dr(K + K — crS m i n ) One can now introduce the Hamiltonian form of 
the path integral via 

j Dxe^ CdT = J DxDpe ipk±k -^ HdT 
and rewrite the first two exponents as 



exp ji J PiRidr + ij[-p i (l- 2x(t)) + q^dr 



The term proportional disappears because of boundary conditions r M (0) = r M (T) = 0, and 
one obtains 

p + = l -P + {l-2x), [p + ,r_] = -t 
One can rewrite this in the form 

d . 

Separating out the center of mass motion one defines the Hamiltonian from the correspond- 
ing Minkowskian action A M : 



J dz + L M , H = p ± r ± - L M (165) 
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with p± = dL M /df±. 

After all one easily obtains the explicit form of the Hamiltonian 

1 (ml ml L 2 (j )±r± j rir _y 

H = o i 1 1 2 ^ ~ 

J v Hi + r\ J 

Let us remind, that /ii, /i 2 are einbein fields, playing the role of P + momenta of the particles, 
jl = 111/22/(^1 + ^2) , ^0 = Jo d(3is(f3) and v((3) is the einbein field with the physical meaning 
of the P + -momentum, carried by the string. As it has been mentioned above, the Hamiltonian 
explicitly depends on L 2 = p \r \ — (p±r±) 2 , the corresponding mass parameter is equal to 

(HI): 

1 

\ 2 1 2 1 f JQ,,(n\(n \2 



a = Hi(l — x) + jj, 2 x + / dj3v((3)(l3 — x 





The variable x in the above expression is the same as defined in ( |161| ): 

x = (167) 

where 

P+ = in + fi2 + u (168) 
is the light-cone total momentum of the system. Also 

7 = vo ((/3) - -^-) (169) 

where ((3) = Jq 1 df3f3v(f3)/ Jq 1 df3v(f3). It is easy to show, that in the limit a <^ m 2 Hamiltonian 
turns out to have well known form: 

1 /m? ml \ , , , 

H =2?vU 1 + T^) +ff - |r ' (170) 

Let us stress that separation of this kind cannot be done beyond nonrelativistic limit. 

As the next step one should quantize the classical Hamiltonian function. Before doing it, 
one should choose the appropriate set of dynamical variables. Three einbein fields Hi,fi2,v 
introduced above play different dynamical roles. In the nonrelativistic case m\,m2 ^> \fo (and 
therefore for the free particles) the dependence of Hamiltonian (and wave functions) on v can 
be correctly found by minimization procedure with v taken as a classical variable. This is in 
its turn a consequence of the fact, that string in our approach is the minimal one - it has no 
internal degrees of freedom and may only stretch or rotate as a whole. 

On another hand, fi\ and H2 on the light cone play the role of legitimate quantum dynamical 
degrees of freedom and can be expressed through x and P + as in ( |168| ) and ( |167|) . 

There are two canonically conjugated pairs in our problem: {p±,f±} and {x, (P + r_)}. We 
introduce also a new dimensionless variable y instead of v : y((3) = v(f3)/P + . It satisfies the 
obvious condition: < y < 1. This variable depends on f± as well as on (P + r_). Rigorously 
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speaking one should extract this dependence by the minimization of the Hamiltonian with 
respect to y as it has been explained above. Instead an easier (however approximate) way is 



chosen. On physical grounds it can be shown that for small rj_ one has yo = J ydf3 ~ const ■ rj_ 
so that linear string energy density yo/rj_ stays finite if r\ — » 0, and using this property one 
can reproduce the correct 1+1 limit, namely 't Hooft equation [72], from the 3+1 light-cone 
Hamiltonian ( |166| ) (see ffl |. On the other hand, if distances are increasing, y~o tends to some 
limiting value which is determined by the virial theorem arguments. So one can parameterize 
y introducing several parameters and replace the minimization in the functional sense by the 
ordinary minimization with respect to these parameters. We have chosen the simplest 2- 
parametric form for y: 

V=^~ (m) 

where t = r\ and y and a are free parameters. The requirement < y < 1 leads to the 
restriction < y < a. Let us stress again that this parameterization is the matter of convenience 
and all physical results are determined from the requirement that every energy level should have 
its own minimum. 

It is easy to see from (|167|) , that the x-variable is the part of the total momentum, carried 
by the first quark itself and a part of the string, "belonging" to this quark. Rewriting ( |167|) in 
the form: 

1 - x) - y(l - </3» (172) 



x 



v(B) ; f 



one can conclude, that for the given y the variable x may vary in the following limits: 

y(P) <x<l-y(l- ((3)) (173) 
It is more convenient to make rescaling to a new variable p which varies from zero to unity: 



x = y((3) + (1 - y)p ■ < p < 1 
The quantity we are really interested in is the mass operator squared: 

M 2 = 2 H 



(174) 



(175) 



Using definition ( |167| ) we obtain the following form for M 



2. 



M 



+ 



1717, 



+ " 



+ 



x-Vo(fl) ' l-a;-j/o(l- (/?)) ' cr\ 

1 1 



x - yo((3) 



x-y (l-(P)) 



\p±r± + l{P+r- 



here 



+ / C r\d(5 + - 
J v 1 



Vo (P+r- 



y i - yo r\ 

c=(x-y ((3))a-x) 2 + 
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[l-x-yo(l-(P))W+ / y((3-x) 2 dp 

Jo 

and 

?=T^«0-*) 

As it was mentioned above, it is more convenient to express M 2 through a new variable p: 

1 - Vo 

In terms of this variable one obtains: 

» r ? 1 f m? ml L 2 7 
1 - l/o 1 P 1 - P c t 

\p 1-p) t 

r 2 



(ire) 



where the notation £ = r[ was used. 

Quantization of the above expression is done according to the canonical commutation rela- 
tions: 

{p k ± ,r{} = -i6 kj 
{x, (P+r-)} = -i 

We are looking for the wave function of the problem given in the mixed coordinate - momentum 
representation if) = ip(p,t), so one has to substitute into ( [L76|) the operators: 

(P +r _) = i (—^—) ; P± = -i ° 



1 — yo J dp ' dr\ 
The important point is the operator ordering. We use the Weil ordering rule, i.e. 

AB -> ^{Ab + bA) 

for any noncommuting operators A and B. Let us notice that y explicitly depends on t according 
to ( |171| ) and hence should also be differentiated during the ordering procedure. The final result 
for the operator M 2 may be found by the straightforward calculations, it is 

,-.,2 . o 2 , s 2 , a 2 

+A 4t +A 4 P +A ° (i77) 



where A t are expressed through p,t,y, (see for more details). 

We have chosen six sets of quark mass parameters, including 4 sets of equal masses and 
2 sets of unequal masses in the interval from to 5 GeV. The values of parameters y, a, e 
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are obtained by minimization procedure as explained above. The eigenstates obtained by our 
numerical procedure are shown in Fig.2. 

Let us first discuss the light-cone spectrum. The Hamiltonian ( |166| ) depends explicitly on 
L z , whereas the Lorentz-invariance requires that masses do not depend on L z , but depend on 
L. Consequently one expects that the lowest state for L z = corresponds to L — 0, while the 
next state with L z = corresponds to L — 1 and must be degenerate in mass with the state 
L z = 1, L — 1. This is clearly seen in the data. 

At this point one should take into account another degeneration - the dynamical one, and 



to this end discuss first center-of-mass spectrum, obtained from the cm. Hamiltonian [37 



The cm. Hamiltonian for L = reduces to the so-called spinless Salpeter equation [[71] which 
was actually solved by us. The approximate form of the spectrum can be represented as: 

L, N r )Y = 2ncr ^2N r + ^-L + ^\ (178) 

For L > one should take into account the string contribution, absent in spinless Salpeter 



equation but present in the Hamiltonian |)7|]. This correction is done in (113). The corrected 
masses, M(L, N r ) = M^°> + AM are shown in Fig.2. Now one can see that there is an approxi- 
mate degeneration in mass M = M^°> + AM of states when one replaces one unit of N r by two 
units of L. The same type of approximate mass degeneration is seen in the light-cone spectrum 

- compare e.g. the states with (N r , L, L z ) = (1, 0, 0) ; M LC = 2.25 and (0, 2, 0) ; M LC = 2.28. 

In Fig.2 this degeneration is visualized as the fact that masses appear on the vertical lines. 
This degeneration is a dynamical one and is a characteristic feature of nonrelativistic oscilla- 
tor. In our relativistic case it reveals a new string-like symmetry, typical for the QCD string 
spectrum [El]. 

Let us now compare the light-cone and the center-of-mass spectrum. One could expect the 
coincidence up to an overall shift due to the different treatment of Z-graphs in two systems. 
In the cm. Hamiltonian these Z-graphs are presented but supposed to be unimportant on the 
grounds, that the backtracking of a quark in time necessarily brings about a folding in the string 
world sheet which costs a large amount of action and is therefore suppressed. The situation is 
different in the light cone - it is general belief that Z-graphs are absent here. One also expects 
that the Z-graphs (and the overall shift) should decrease if quark masses are increasing. 

The comparison of the spectra can be made from Fig. 2. One can see indeed some overall 
shift down in the cm. spectrum by some 0.1 GeV and otherwise the masses coincide within the 
accuracy of computation. This fact is highly nontrivial since two quantum Hamiltonians (the 
light-cone and center-of-mass ones) are different, they cannot be obtained from each other by a 
simple boost or other simple transformation. The light-cone Hamiltonian is rather complicated 
and it took the authors more than a year to get reasonable numerical results for it. 

We have also checked the quark mass dependence of the overall shift of spectra and proved 
that it drops sharply with the quark mass increasing supporting the idea, that the shift is due 
to different treatment of Z graphs (or self energy graphs). This fact also confirms, that the 
shift is not a consequence of some systematic errors of our procedure. 

We now turn to eigenf unctions. One expects in this case two types of excitation: the radial 
one leading to new nodes on p coordinate and similar to the 1+1 excited states and the r±- 
excitation which causes nodes in the t-coordinate and associated with the orbital excitation. 
This is clearly seen in Fig. 3, where (a) refers to the ground state, (b) - to the orbital and (c) 

- to the radial excitation. A more complicated example, combining both types of excitation is 
given in Fig. 3 (d). 
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As the next illustration we show in Fig. 4 the case of two heavy quarks, demonstrating two 
types of excitation and also a new feature - the actual region of parameters is squeezed to a 
small region near p = 1/2. 

At this point it is important to find connection of our light-cone wave function to the 
nonrelativistic one, usually defined in the cm.. In [5] it was demonstrated that this connection 
can be established only when both quarks are heavy, mi,m 2 y/a. In this case y <C 1 (as 
can be found by direct minimization of the Hamiltonian (|166|) ) - the string transforms into 
the potential and loses its material and momentum contents. One can introduce as in [5] the 
relative momentum p z and relative coordinate r z : 

p z = (mi + m 2 ) ( x — ) (179) 

V nii+nioJ 



nil + ^2 
nil + ^2 



P r 

r z = ± — (180) 



and the M operator can be written as: 



M « mi + m 2 + — (p I + p I) + a Jr\ + r\ (181) 
2m v 

Hence one can write the momentum-space nonrelativistic wave function ty(p 2 ) directly 
through the light-cone variables: 



#(p 2 ) = $ 



■> \2 / ^1 

p ± + (nil + ™2) x 

nil + m 2 



(182) 



This representation is valid in the large mass limit m« ^> y/a stated above and in addition 
near the center of the x-distribution, i.e. when \x — mi/(mi + m 2 )| 1. The width of the 
peak in x variable is proportional to (m% + m 2 )~ 2 and is very narrow for heavy quarks. The 
form (|179|) is not correct for x at the ends of the interval, i.e. x = 0, 1 (remember that for large 
fflj the extremum value of y tends to zero and x = p changes in the interval [0,1]). Indeed the 
exact wave function as discussed in Section 4 vanishes linearly at x — 0, 1, while the r.h.s of 
( |182j ) stays nonzero. Moreover, the Jacobean J of the phase space d 3 p = Jdxd 2 p± is constant 
J = nil + ^2 and does not change this conclusion. 

The correspondence of the cm. and light-cone wave functions is lost if quark masses 
are of the order of \fo or less. The physical reason is that the role of dynamics is now 100% 
important and the dynamics is different in different frames: the light-cone Hamiltonian and 
wave functions are connected with those of the center-of-mass by a dynamical transformation, 
which includes nonkinematical Poisson operators. Hence one cannot hope to obtain one wave 
function from another using only transformation including the kinematical (free) part of the 
operator M 2 . Moreover, a glance at the operator M 2 in (|166|) helps to realize that the separation 
of M 2 into the purely kinetic part (containing only momentum and masses) and purely potential 
one (containing only relative coordinates) is impossible at all: e.g. the string enters through 
the term (P + r_) 2 and mixed terms like p±r± are present everywhere. This circumstance limits 



the use of simple recipes known in literature [73], which connect the cm. and light-cone wave 



functions. In particular, the ansatz for p z suggested in [|73| and used for equal quark masses 



Pz 



"* + # (183) 



\x(l-x)\ 2 
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coincides with our form ( |180| ) rigorously derived in [|(J only for large m, m ^> -y/cr and x in 
the narrow region around x = 1/2 but differs at the tails of the wave function. For light quark 
masses, m < a and small p± the relation ( |183| ) yields incorrect results which can be seen in 



the almost constant behaviour in the interval [0,1] of the light-cone wave function, produced 
by the insertion of (|183|) into the cm. wave function, typically 



tp(p 2 ) ~ exp(-aV 2 ) -> exp(-a 2 [pi + ^^^ (x- i) 2 ]) 

is insensitive to x for m,p± — > 0, whereas the exact light-cone wave function is decreasing as 
x(l — x), see Fig. 5. 

At the same time for heavy quark masses, m ^> y/a the two functions, one transformed by 
( |183j ) from the cm. and another is a genuine light-cone wave-function solution of (177), are 
very close to each other. 

We now turn to the formfactor of computed states. The main feature is a very slow decrease 



of F(q 2 ) with q 2 which signals in particular a small radius of states. Indeed the values of J (r 2 ) 
are too low (~ 0.338 fm for massless quarks). This fact is in agreement with the earlier cm. 



calculations of [74 



Finally we turn to the quark-distribution function q(p). It is computed through the light- 
cone wave-function ip(p, t) using the relation q(p) = n / °° |\I/(p, t)\ 2 dt and shown in Fig. 6. One 
can see the symmetric behaviour of q(p) with respect to reflection p — > 1 — p. At the ends of 
the interval q(p) vanishes like p 2 and (1 — p) 2 , in the agreement with the 1/q 3 behaviour of the 



formfactor at large q due to the Drell-Yan-West relations [75 . 

Note the narrowing of the peak in q(p) in Fig. 6 for increasing quark masses. 

The main physical idea of our approach is that the most part of nonperturbative dynamics 
in QCD is due to the QCD string, and the latter is described by the Nambu-Goto part of the 
Hamiltonian, which was written before in the cm. [5] as well as in the light-cone coordinates 

H- 

Only valence part of the Fock's column was considered above in the paper, also for simplicity 
spins and perturbative gluon exchanges are neglected. To do the systematic comparison with 
experiment all these three simplifications should be eliminated. Let us discuss their effect 
point by point. The higher Fock states are necessary to reproduce the Regge behaviour of 
q(x) ~ £- a p(°) at small x (and at x — > 1 for quark distributions of hadrons in high-energy 
scattering). Here comes the first crucial point; to be answered in the second part of these 
lectures, what is the QCD reggeon? 

In our method the higher Fock states, constituting the QCD reggeon, correspond to sev- 
eral gluons propagating in the nonperturbative background and therefore confined to the ex- 
cited Nambu-Goto surface [2]. These states are in one-to-one correspondence with the excited 
Nambu-Goto string states. This is the picture at large distances; at small distances smaller 
than the vacuum correlation length (width of the Nambu-Goto string) T g ~ 0.2 fm, the string 
disappears and the usual perturbative gluon exchanges reappear. 

The effect of spin of light quarks is highly nontrivial [3] . It leads to the creation of the new 
vertex, which yields the constituent quark structure. Physically one may imagine this structure 
as being due to the light-quark walks around the end of the string. 

Having said all this, what are the lessons of the present work and of it possible development? 

The first lesson is that valence quark component can be successfully dynamically computed 
on the light cone; the minimal string on the light cone is physically and mathematically well 
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defined. The structure of the spectrum obtained on the light cone for the first time reasonably 
coincides with that of the cm. Hamiltonian for the string with quarks. Effects of Z-graphs 
are estimated to be less than 12% even for light quarks. Moreover, the nonperturbative wave 
function obtained directly on the light cone allows to calculate nonperturbative contributions 
to the formfactor and structure function. 

The second lesson is that the formfactor computed directly on the light cone is close to 
of cm. for small q, but is systematically above the cm. formfactor for larger q. This is not 
surprising, since in the light-cone formfactor there is a mechanism of the "redistribution" of 
the momentum q between the quarks, since q enters the light-cone formfactor multiplied with 
(1 — x), so that the larger q, the smaller is (1 — x) and the wave function does not decrease too 
fast. Physically it means that at large q the configuration survives where the spectator quark 
gets as little momentum as possible so that it can be easily turned together with the active 
quark. This is exactly what is called the Feynman mechanism [76]. 

The third lesson is that the minimal string plays only a passive role on the light cone, 
namely it participates in sharing of the total momentum and carries the part equal to (y) , but 
it does not produce the x-distribution in structure function, which could simulate the gluonic 
structure function. The reason is that the string variable - the einbein field v - is quasiclassical 
and has no dispersion. 

The value of {y) computed according to relation < > = n J °° Jq 1 \^/(p,t)\ 2 dpdt depends on 
quark masses and is equal to 0.22 for massless quarks. It is resonable that (y) = 0.22 is smaller 
than the experimental value of overall gluon momentum, 0.55, since in our picture the difference 
should be filled in by higher Fock components. 

The fourth lesson comes from the comparison of the computed quark distribution, with the 
experimental data for the pionic structure function JT7|. Behaviour of q(p) at small p and small 
(1 — p) is symmetric in Fig.6, while in reality q(x) should rise at small x like x~ ap ^ ~ x~ 5 
(we neglect at this point the difference between x and p, which is due to (y)). This peak 
at the small x should be filled in by the contribution of higher Fock components, containing 
additional gluons on the string, as was discussed above. The behaviour of q(p) at p — 1, which 
is calculated to be (1 — p) 2 will be also changed into (1 — p) due to gluon exchanges, which 
account for the formfactor asymptotic 1/q 2 at large q, and the Drell-Yan-West duality ensure 
the (1 — p) behaviour around p = 1. 

Finally, the formfactor calculated above in the paper, shows too little radius of the "pion" 
(r 2 ) pa (0.34 fm) 2 as compared with the experimental one (r 2 ) ~ (0.67 fm) 2 . This fact is in 



qualitative agreement with other calculations, where the cm. wave-function was used [74|, and 



some authors assumed as in |7(| that quarks have "internal" structure and their own radius 
which should be added to the "body radius" to reproduce the experimental value. This fact of 
small body radius seems to be a necessary consequence of the simple string + point-like quarks 
picture, and probably cannot be cured by the higher Fock components;in the second part of 
lectures we shall come to this point fromanother direction. 



8 Conclusions 

As one can see from previous chapters, the simple relativistic string Hamiltonian (111-112) and 
its extension to the hybrid and glueball case works surprisingly well and describes quantitatively 
the observed meson spectrum and calculated on the lattice spectrum of hybrids and glueballs. 
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This implies that confinement in the form of the thin string is the basic element of the strong 
interaction, and for the most hadrons this is enough to reproduce main features like masses, 
decay constants, magnetic moments etc. Last years there appeared new successful applications 
of the method to the polarizability of mesons [78], to the description of nucleons [79] and to 
the calculation of the nonperturbative correlation masses in SU(2)- Higgs model [80]. Analysis 
of heavy quarkonia, sensitive to small distances adds to this picture two elements: one needs 
to know one-loop corrections to spin-dependent forces and small distance behaviour of NP 
forces [19], [30]. However the world of Nambu-Goldstone bosons and chiral symmetry breaking 
(CSB) phenomena was completely omitted in the previous material for the lack of space. A 
new approach based on the analysis of the heavy light system [3] , allows to formulate dynamical 
equations for mesons and baryons [65] in this case, taking into account CSB. Simultaneously 
one can calculate in the method the NP self-energy of the quark (constant Cq subtracted from 
the mass eigenvalue in (112)), and consider spin effects not assuming its smallness. This part 
of material will be published separately as a second part of lectures. 

The author is grateful to A.B.Kaidalov, Yu.S.Kalashnikova, F.J.Yndurain for discussions 
and to the Organizing Committee and Professor Lidia Ferreira for the effective organization 
and warm hospitality at the XVII International School. 
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Table 1 



Comparison of the WKB spectrum of Hamiltonian (72) with the exact spectrum of 

Hamiltonian (71) for m = and a = 0.2 GeV 2 . 



n 





1 


2 


3 


4 


5 


n 


1.475 


2.254 


2.825 


3.299 


3.713 


4.085 


n 


1.412 


2.106 


2.634 


3.073 


3.457 


3.803 



Table 2 

Effective mass eigenvalues Ho(n, I) (in GeV for 07 = 0.18GeV 2 ) obtained from Eq.( 70), (72) 
A*o = y/^Ji^f^) 3 ^- upper entry,and eigenvalues of reduced equation a(n)-lower entry. 



n 





1 


2 


3 


L 













0.352 


0.535 


0.67 


0.78 




2.3381 


4.0879 


5.520 


6.786 


1 


0.462 


0.611 


0.732 






3.3613 


4.8845 


6.216 




2 


0.55 


0.68 








4.248 


5.63 






3 


0.627 










5.053 









Table 3 

Mass eigenvalues of the rotating string Hamiltonian (67) 
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n I 


1 


2 


3 


4 


5 





1.865 


2.200 


2.481 


2.729 


2.956 


1 


2.562 


2.832 


3.068 


3.281 


3.480 


2 


3.091 


3.329 


3.540 


3.733 


3.913 


3 


3.535 


3.753 


3.947 


4.125 


4.290 


4 


3.925 


4.128 


4.309 


4.476 


4.629 


5 


4.278 


4.469 


4.638 


4.797 


4.939 



Table 4 

Energy eigenvalue ct(A) of the reduced equation (87) as a function of the dimensionless 
paramenter A, (88); derivative |a'(A)| =| | and wave function at origin xa(0) (solution of 
(87)) divided by xo(0) (solution at A = 0) 



A 





0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1 


a(A) 


2.338 


1.99 


1.896 


1.801 


1.704 


1.604 


1.502 


1.398 


|a'(A)| 


0.84 


0.9 


0.92 


0.94 


0.96 


0.98 


1.01 


1.04 


Xa(0) 
Xo(0) 


1 


1.249 


1.318 


1.392 


1.469 


1.549 


1.633 


1.72 



Table 5 

The factor p(AF), as a function of paramenter A, (88), for two values of Aq C d entering (88) 
via one-loop expression (17) 

_ Xaf(0) 
P X(0) 



Aqc*d;A 





0.4 


0.6 


0.7 


0.8 


1 


90 MeV 


1 


0.976 


0.92 


0.884 


0.838 


0.76 


140 MeV 


1 


0.984 


0.95 


0.911 


0.877 


0.79 



Table 6 

Mass correction (—AM) in GeV, Eq. (112), due to the rotating string, computed quasiclas- 
sically. The entries in parenthesis are obtained from Eq. (113) 
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n 





1 


2 


L 








1 


0.082 
(0.066) 


0.014 


0.005 


2 


0.119 
(0.120) 


0.035 




3 


0.164 
(0.167) 
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Table 7 

Meson masses for the Hamiltonian (85) with the parameters 

a s = 0.39 

a = O.UGeV 2 
C = -482MeV 



meson 


mass (theory) in MeV 


mass (experiment) in MeV 


p(lS) 


769 (fit) 


768.5 ±0.6 


p(2S) 


1500 


1465 ± 25 


p(3S) 


2035 


2149 ± 17 


n(2S) 


1325 


1300 ± 100 


tt(35) 
M(l 3 Pj) 


1895 


1861 ± 13 


1200 


1230 ±22 = M cog (i 3 P,,) 
candidates") 


M'2( 3 Pj) 


1790 


a 2 (1700),a 1 (1700) 


M"3( 3 P,) 


2265 


ai(2100),a (2050) 



a ) T.Barnes |hep-ph/9907259| 



Table 8 



Compilation of lattice results for exotic 1 
tions of the present approach 



hybrid masses (in GeV) compared to calcula- 





M(bbg) - M(bb) 


M(ccg) - M(cc) 


M(ssg) 


M(udg) 




1.14[a] 


1.34 [e] 


2.00[g] 


1.88[g] 




1.3[b] 


1.22 [f] 


2.17[e] 


1.97[e] 


Lattice*) 


1.54[c] 


1.323[c] 




1.90[i] 




1.49[d] 






2.11[f] 


Present 
approach 


1.5 


1.45 


2.34 


2.16 



) From D.Toussaint, |hep-ph/9909088| ; [a],...[g] refer to sources in the paper. 
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Table 9 



Lattice data on glueball masses (in GeV) compared to the present approach [60] 





[61] 


[62] 


[63] 


Present approach [60] (07 = 0.23GeV 2 ) 


0++ 


1.73 ±0.13 


1.74 


1.69 


1.58 




2.67 ±0.31 


3.15 


2.48 


2.71 


2++ 


2.4 ±0.15 


2.48 




2.59 






3.22 




3.73 


o-+ 


2.59 ±0.17 
3.64 ±0.24 


2.38 


2.54 


2.56 
3.77 


2-+ 


3.1 ±0.18 
3.89 ±0.23 


3.38 


3.31 


3.03 
4.15 


3++ 


3.69 ±0.22 


4.31 


4.27 


3.58 


1— 


3.85 ±0.24 






3.49 


2— 


3.93 ±0.23 






3.71 


3— 


4.13 ±0.29 






4.03 


1+- 


2.94 ±0.17 


3.03 






2+" 


4.14 ±0.25 




4 




o+- 


4.74 ±0.3 









Table 10 

Characteristics of heavy-light mesons in comparison with experiment and lattice data 



meson 


B 


B* 


B s 


D 


D* 


D s 


D* 


M-m Q 


0.479 


0.539 


0.535 


0.496 


0.630 


0.532 


0.666 


M 


5.279 


5.339 


5.335 


1.896 


2.030 


1.932 


2.066 




5.279 


5.324 


5369 


1.869 


2.010 


1.968 


2.112 


f M (Gev); 
= 0.427 


0.183 






0.221 




0.264 




Mat) 
JM 

= 0.427 


0.210 




0.251 
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Figure 1: Quasiclassical spectrum of Hamiltonian (67) for m = and a = 0.2GeV 2 . The 
leading experimental Regge trajectory in angular momentum L is given in the upper plot for 
comparison. Theoretical prediction for the p-meson mass, with color Coulomb and spin effects 
included is shown by the dot at L = and does not violate the straight-line behaviour of the 
leading trajectory 



Figure 2: The Chew-Frautschi plot with masses computed via the light-cone (circles, squares 
and triangle for L z = 0,1,2 respectively) and the cm. Hamiltonian (stars). The systematic 
overall mass shift is seen as a divergence of straight lines passing through circles and stars. 
The states with high L or high N r (daughter trajectories) are numerically less reliable and not 
shown. 



Figure 3: The 3d plots of wave functions of the four lowest states of light-cone Hamiltonian 
for zero quark masses. Coordinates on horizontal plane are < p < 1 , < t < 15 (in units of 



Figure 4: The same as in Fig. 3 but for heavy quark masses, mi = mi = 5 GeV. 



Figure 5: The 3d plots of the ground-state wave functions \l/(p, t), computed via the light-cone 
Hamiltonian (upper part) and via the cm. Hamiltonian, with the standard substitution (|183| ) 
(lower part) for zero quark masses. 



Figure 6: The quark-distribution function q(p) computed with light-cone wave-functions for 
the cases considered. 
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